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Abstract 

Neutron stars provide a fertile environment for exploring superfluidity under 
extreme conditions. It is not surprising that Cooper pairing occurs in dense 
matter since nucleon pairing is observed in nuclei as energy differences be- 
tween even-even and odd-even nuclei. Since superfluids and superconductors 
in neutron stars profoundly affect neutrino emissivities and specific heats, their 
presence can be observed in the thermal evolution of neutron stars. An ever- 
growing number of cooling neutron stars, now amounting to 13 thermal sources, 
and several additional objects from which upper limits to temperatures can be 
ascertained, can now be used to discriminate among theoretical scenarios and 
even to dramatically restrict properties of nucleon pairing at high densities. In 
addition, observations of pulsars, including their spin-downs and glitch histories, 
additionally support the conjecture that superfluidity and superconductivity are 
ubiquitous within, and important to our understanding of, neutron stars. 



1 Introduction 

In this contribution, we describe the roles of neutron superfluidity and proton 
superconductivity in the astrophysical setting of neutron stars, drawing upon 
lessons learned from similar phenomena occuring in laboratory nuclei. We will 
focus on both the thermal evolution (i.e., cooling) as well as the dynamical 
evolution (i.e., spin-down) of neutron stars. In the former, pairing dramatically 
affects neutrino emission processes and the specific heat of dense matter. In 
the latter, pairing may be responsible for the observed anomalous values of the 
braking index and the glitch phenomenon. We will also briefly describe the 
possibility of pairing in the presence of hyperons and deconfined quarks. 

In fermionic systems, superfluidity and superconductivity occur due to the 
pairing of neutral and charged fermions, respectively. The Cooper Theorem [1] 



states that, in a system of degenerate fermions the Fermi surface is unstable 
due to the formation of "pairs" if there is an attractive interaction in some spin- 
angular momentum channel between the two interacting particles. The essence 
of the BCS theory [5] is that as a result of this instability there is a collective 
reorganization of particles at energies around the Fermi energy and the appear- 
ance of a gap in the quasi-particle spectrum. This reorganization manifests itself 
in the formation of "Cooper pairs" . At high-enough temperature, the energy 
gap disappears and the system reverts to its normal state. 

To begin, in Sec. [2j we motivate the existence of pairing in neutron stars by 
examining the pairing phenomenon in laboratory nuclei. We then summarize 
the relevant properties of neutron stars, including their interior compositions 
and properties of their crusts in Sec. |3l We describe in Sec. |4] how pairing in 
dense matter is achieved and, in Sec. [5l we present a brief description of ex- 
pectations of pairing in deconfined quark matter that may be present in the 
inner core of the most massive neutron stars. Initially, neutron stars cool pri- 
marily due to neutrino emission in their interiors before surface photon cooling 
takes over later in their lives. In Sec. [6j we summarize the various neutrino 
processes that can occur in dense matter. Theoretically, these processes can 
proceed either very rapidly (enhanced neutrino emission) or relatively slowly. 
We also describe an important secondary process that greatly influences the in- 
terpretation of observations - neutrino pair emission from the pair breaking and 
formation (PBF) of Cooper pairs. This process is triggered as the ambient tem- 
perature, decreasing because of cooling, approaches the critical temperatures for 
superfluidity and superconductivity. The discussion in this section shows how 
the existence of superfluidity or superconductivity dramatically influences neu- 
trino emissions, leading to both the quenching of enhanced neutrino emission 
and bursts of neutrino emission due to the PBF processes. 

We present in Sec. [7] simplified analytical models of neutron star cooling in 
order to gain physical insight. These analytical models are complemented by de- 
tailed numerical simulations which include general relativity and state-of-the-art 
microphysics, such as the dense matter equation of state, thermal conductivi- 
ties, neutrino emissivities, and specific heats. We also summarize in this section 
the abundant observational data, consisting of estimates of surface temperatures 
and ages, which collectively describe the thermal evolution of neutron stars. We 
will show in Sec. |8] that the bulk of this data supports the so-called "Minimal 
Cooling Paradigm" , which supposes that no drastically enhanced neutrino emis- 
sion processes occur, or if they do, they are quickly quenched by superfluidity 
or superconductivity. Nevertheless, a few sources suggesting enhanced cooling 
are observed, and we discuss the implications. 

In Sec. [9] we describe an outstanding recent development in which the first 
real-time cooling of any isolated neutron star - the young neutron star in the 
supernova remnant of Cassiopeia A - is observed. The observed rate of cooling 
is more than 10 times faster than expected, unless both neutron superfluidity 
and proton superconductivity are present in the star's core. These observations 
provide the first direct evidence for superfluidity or superconductivity in the 
interior of a neutron star, and can be verified by continued observations of the 



2 



neutron star in Cassiopeia A. 

Section 10 is devoted to other observations of neutron stars and their dynam- 
ical evolution that may also indicate the presence of superfluidity. One concerns 
the deceleration observed in the spin-down of pulsars, which could be due to 
superfluidity in their cores, while the other is related to sporadic spin jumps, 
commonly known as glitches, thought to stem from superfluidity in neutron star 
crusts. 

Summarizing discussions and conclusions are contained in Sec. 11 



2 Pairing in Nuclei 

Soon after the development of the BCS theory, Bohr, Mottelson & Pines [3] 
pointed out that excitation energies of nuclei exhibit a gap, as shown in the 
left panel of Fig. [T] A nucleon in even-even nuclei, whether neutron or proton, 
clearly requires a minimum energy for excitation. This energy was interpreted 
as being the binding energy of a Cooper pair which must break to produce an 
excitation. In contrast, odd-even nuclei do not show such a gap, and this is due 
to the fact that they have one unpaired nucleon which can be easily excited. 
The right panel of Fig. [l] shows that pairing also manifests itself in the binding 
energies of nuclei, even-even nuclei being slightly more bound than odd-even or 
odd-odd nucleiil 
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Figure 1: Left panel: Lowest excitation levels of nuclei (adapted from ^J). Right 
panel: Binding energy per nucleon for the most beta-stable isobars (adapted 
from [4]). 

The systeniatics of neutron pairing energies in nuclei, defined through 

^Z,N = — '2Ez,N + Ez,N-l) , (1) 

^Notice that, as a result of pairing, the only stable odd-odd nuclei are ^}1(Z = 1,N = 1), 
^Li(3,3), ^"3(5,5), and ^*N(7,7). All heavier odd-odd nuclei are beta-unstable and decay into 
an even-even nuclei. 
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Neutron number N 

Figure 2: Upper panel: Odd-N pairing energies. Lower panel: Even-N pairing 
energies. Figure adapted from [5]. 

where Ez,n is the binding energy for charge Z and neutron number N, and 
+(— ) refers to odd-N and even-N nuclei, are shown in Fig. |2] A few key facts 
to note are: 

1. Pairing energies range from about 3 to 0.5 MeV, decreasing in magnitude 
with increasing neutron numbers; their behavior with the mass number 
A^N + Z is well fit by 

An,z = 24/^ + 0.82 ±0.27, for TV odd, (2) 
^N,z = 41/A + 0.94 ±0.31, for TV even. (3) 

2. Dips (peaks) occur adjacent to (at) the sequence of magic numbers 
N = 14, 28, 50, 82 and 126 for iV odd (even). 

Systematics of proton pairing energies for odd- and even-Z as a function of 
proton number Z (see Fig. 2 and Table 1 of [5]) show similar qualitative be- 
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havior, but the magic number effects are less pronounced. If the pairing energy 
were to be extrapolated to infinite matter using Eq. ([2| and A ^ oo, the gap 
should vanish. As gaps in infinite matter are predicted to be finite, a saturation 
penomenon is at play. Extended matter in '^He, albeit for different reasons, also 
exhibits finite gaps. 

In addition to the excitation spectra and binding energies of nuclei, the pair- 
ing phenomenon plays important roles in the dynamical properties of nuclei such 
as their rotational inertia and the large amplitude collective motion encountered 
in fissioning nuclei. Tunnelling effects in spontaneously fissioning nuclei receive 
an enhancement factor 2A^/G'^, amounting to an order of magnitude or more (G 
is a typical pairing interaction matrix element between neighboring mean-field 
configurations) . 

Pairing effects are also evident in nuclear reactions. For example, thermal 
neutrons of energy only ~ 0.025 eV are needed to cause the fission of ^g^U 
(which results in the even-even compound nucleus ^12^), whereas fast neutrons 
of higher energy ~ 1 MeV are needed to induce fission of ^|fU (the compound 
nucleus in this case is the even-odd nucleus ^gfU) . What is interesting is that this 
phenomenon was appreciated well before the BCS theory was formulated, and 
it lies at the root of building nuclear reactors and purifying naturally-occurring 
uranium to contain more of the 235 isotope than the 238 isotope. 

Besides BCS pairing, the pairing energies shown in Fig. [2] receive contribu- 
tions from other sources since nuclear sizes are much smaller than the coherence 
length of the pairing field. The odd-even staggering is caused by a combina- 
tion of effects including the pair-wise filling of the orbitals, diagonal matrix 
elements of the two-body interaction, three-nucleon interactions, the bunching 
of single particle levels near the Fermi energy, and the softness of nuclei with 
respect to quadrupolar deformations. The global description of the pairing phe- 
nomenon in nuclei is based on the Hartree-Fock-Bogoliubov approximation and 
recent accounts may be found in Refs. [Sj [6l [7]. The basic cause for pairing in 
nuclei is, however, easy to identify. The nuclear interaction between identical 
nucleons is strongly attractive in the spin S = channel, the di-neutron being 
nearly bound. The even stronger attraction between neutrons and protons in 
the spin 5 = 1 channel produces bound deuterons, but its effects are mitigated 
in heavy nuclei due to the imbalance of neutrons and protons and attendant 
many-body effects. In any case, with proton and neutron pairing energies on 
the order of an MeV, nuclei represent the highest temperature superconductors 
and superfluid objects in the laboratory. It is interesting that gaps of similar 
order-of-magnitude are expected for nucleon pairing in neutron stars. 

3 Neutron Stars 

Given that nucleon pairing is important in nuclei, we should expect that pairing 
will also occur within neutron stars, as was originally pointed out by Migdal in 
1959 [5]. Although matter within neutron stars may be heated to more than 
10^^ K during birth, and may remain warmer than 10^ K for hundreds of thou- 
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sands of years, the nucleons are generally extremely degenerate. Furthermore, 
given the high ambient densities, the critical temperatures for pairing to occur 
are large, ^ 10* — 10^° K. The onset of pairing is expected to take place in some 
parts of a neutron star's interior within minutes to thousands of years after 
birth, and is expected to lead to alteration of several important properties of 
matter. While pairing will not affect the pressure-density relation significantljj^ 
and, therefore, the overall structure of neutron stars, the specific heat of dense 
matter and the emissivity of neutrinos are dramatically influenced. Both emis- 
sivites and specific heats are altered at and below the critical temperature, and 
when the temperature falls well below the pairing critical temperature, both 
vanish exponentially. This has important consequences for the thermal evolu- 
tion of neutron stars that will be described in several subsequent sections of this 
chapter. 

Superfluidity can also important in the dynamical evolution of neutron stars. 
It has long been suspected that the so-called "glitch" phenomenon observed in 
pulsars is due to the existence of superfluids within neutron star crusts and 
perhaps their outer cores. Superfluidity within neutron stars might also signif- 
icantly contribute to the anomalous braking indices, which are related to the 
observed long-term deceleration of the spin-down of pulsars. These phenomena 
will be discussed in Sec. [TOl 

Neutron stars contain the densest form of cold matter observable in the 
Universe, in excess of several times the central densities of nuclei (which is often 
referred to as the nuclear saturation density pnuc — 2.7 x 10^^ g cm~'^). Note 
that pnuc corresponds to the density where cold matter with a proton fraction 

= 1/2 has zero pressure. While larger mass-energy densities are transiently 
reached in relativistic heavy ion collisions, the resulting matter is extremely 
"hot". Black holes contain a much denser form of matter, but their interiors are 
not observable. It has long been believed that neutron stars can only form in the 
aftermath of the gravitational collapse of a massive star [9] , commonly known 
as gravitational-collapse supernovae to distinguish them from thermonuclear 
explosions of white dwarfs leading to Type la supernovae. However, even the 
collapse of stars with masses greater than about 25Mq are thought to produce 
briefly existing proto-neutron stars before they collapse further into black holes 
[TU] . Nevertheless, the vast majority of gravitational-collapse supernovae, due 
to the preponderance of lower-mass stars, will produce stable neutron stars. 

Two simple arguments can convince us that neutron stars, very small and 
very dense, can exist. First, consider the fastest known radio pulsar, Terzan 
5 ad (AKA PSR J1748-2446ad) [11^, and posit that the observed period of its 
pulses, P = 1.39 ms, is its rotational period. (Pulses produced by binaries or 
oscillations of neutron stars are ruled out because nearly all pulsars are observed 
to be slowing down, while orbital and vibrational frequencies increase as energy 
is lost.) Using causality, that is, imposing that the rotation velocity at its 

^At asymptotically-high densities where deconfincd quark matter is thought to exist, pair- 
ing gaps could be of order 100 MeV, in which case the EOS is moderately affected by the 
pairing phenomenon. 
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equator is smaller than the speed of light c, one obtains 



2t:R cP 
v^^ ^riR^ <c ^ R< — = 65 km . (4) 

This value of 65 km for the radius R is but a strict upper limit; detailed theo- 
retical models and observations indicate radii on the order of 12 km. Secondly, 
assuming that the star is bound by gravity, we can require that the gravitational 
acceleration at the equator is larger than the centrifugal acceleration Oeq 
and obtain 

GM ^2 47r2i? M 47r2 

5eq = > acq R= or —> — 

=^ P= 4^ > 8 X 10'' . (5) 

Obviously, Newtonian gravity is not accurate in this case, but we can neverthe- 
less conclude that the central density of these stars is comparable to, or likely 
larger than, Pnuc- Theoretical models show that densities up to lOpnuc [H] 
are possibly reachable. In short, a neutron star is a gigantic, and compressed, 
nucleus the size of a city. 



3.1 The Neutron Star Interior 

A "pure neutron star", as originally conceived by Baade & Zwicky [51 and 
Oppenheimer & Volkoff [T^] cannot really exist. Neutrons in a ball should 
decay into protons through 

n ^ p + +Ve ■ (6) 

This decay is possible for free neutrons as m„ > mp + rrie, where the masses 
denote rest masses. However, given the large densities expected within the 
neutron star interior, the relevant quantities are not the masses, but instead the 
chemical potentials fj,i {i denoting the species) of the participants. The matter is 
degenerate as typical Fermi energies are on the order of 10 — 100 MeV, whereas 
the temperature drops below a few MeV within seconds after the birth of the 
neutron star [14]. Starting with a ball of nearly degenerate neutrons, the decay 
of Eq. ^ will generate a degenerate sea of protons, electrons and anti-neutrinos. 
The interaction mean free paths of anti-neutrinos (and neutrinos) far exceed the 
size of the star. In a neutron star, these can be assumed to immediately vacate 
the star, implying that //^ = 0. Thus, this reaction will ultimately result in the 
/3-equilibrium condition 

^J'7l = + Me , (7) 

which is equivalent to the energy minimization condition de/dxp = 0, where 
e is the energy density and the proton fraction is Xp. At finite, but small 
temperatures, the inverse reaction 

p + e~ ^ n + Ve (8) 
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also occurs since not all particles are in their lowest energy states at all times. 

A neutron star, however, is not born from the collapse of a "ball of neutrons" , 
but rather from the collapse of the iron core of a massive star. At densities typ- 
ical of pre-collapse configurations, p ~ 10^ g cm~'^, Xp > 0.4. During collapse, 
the reaction Eq. ([s]) therefore initially dominates over Eq. (|6| in order to reduce 
the proton fraction. As the density increases, however, the neutrino mean free 
paths become smaller than the collapsing core's size, so Ve temporarily become 
trapped within the core. In this case, /i^^ > 0, altering /3-equilibrium and per- 
mits Xp to remain relatively large. Only after neutrinos are able to diffuse away, 
a time of approximately 10 seconds [14] . will the final /3-equilibrium condition 
Eq. ([t]) be achieved. 

Eq. Q allows one to determine the composition of cold dense nucleonic 
matter. Near pnuc, a rough approximation to the difference of nucleon chemical 
potentials is 

- - 4S'„ Vp7/W (1 - 2xp) , (9) 

where Sy ~ 31 MeV is the bulk nuclear symmetry energy parameter. In dense 
neutron star matter, electrons are relativistic and degenerate, so charge neu- 
trality implies equality between electron and proton number densities and 

fie = he (STT^nsXpY^^ , (10) 

where ub = nn + = p/ms is the baryon density and vriB is the baryon mass. 
Therefore, Eq. ([t]) implies 

yl C \ 3 / \ 1/2 

45^ \ rriB P ^ 



he 



) ^{^) (l-2.pf. (11) 



This is a cubic equation for Xp{p), and at the density pnuc the solution is a;p^nuc — 
0.034. Neutron stars arc therefore composed predominantly, but not completely. 



of neutrons. Eq. (11) predicts that Xp increases roughly with y^p near pnuc- 

Notice that once Pe > ~ 105 MeV, muons will appear and be stably 
present with the conditions of beta equilibrium 

Mm = Me (12) 

and charge neutrality 

Up = rip + He. (13) 

The condition for the appearance of muons is fulfilled when the density is slightly 
above Pnuc- However, even though Hp and are no longer equal, the trend that 
Xp slowly increases with density is not altered by the presence of muons. Fur- 
thermore, in all processes we describe below, there will always be the possibility 
to replace electrons by muons when the density is large enough for them to 
appear. 

At all but the highest densities, nucleons can be regarded as non-relativistic 
in neutron star matter, but muons can be either relativistic or non-relativistic 
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Figure 3: Schematic illustration of the structure of a neutron star; figure taken 
from |15j . The outermost layers of a neutron star, the atmophere, envelope, 
and crust are described in Sec. |3.2[ Superfluidity in the crust is schematically 
represented in inset "A" , and a diagram of the pasta phases in the crust is shown 
in inset "B" . The core is separated into the outer core, which has the structure 
given in inset "C" , and the inner core whose nature is currently unknown. 
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depending on their abundance. For nucleons, the generahzation of the simple 
approximation Eq. Q is 



For muons, 

Mm = ^Jmlc^+Pl^.c\ (15) 

where ppi = fikpi is the Fermi momentum of species i, and Vn and Vp are the 
mean-field energies of n and p. The Fermi momenta are related to the particle 
densities by k\^ = iTr'^m. For the leptons, T4 and are negligibly small. 
With a knowledge of and Vp, the two /3-equilibrium relations Eqs. ([7| and 



( 12 ) can be solved. With four chemical potentials and two equations, a unique 



solution is obtained by imposing charge neutrality, Eq. (12), and fixing ub- 
With the particle densities and chemical potentials known, one can calculate 
any thermodynamic potential, in particular the pressure P and energy density 
e = pc^. Varying the value of ub gives us the equation of state (EOS): P{p). 
Given an EOS, an integration of the Tolman-Oppenheimer-VolkofF (TOV, |13) ^ 
equations of hydrostatic equilibrium provides us with a well defined model of a 
neutron star. 



The potentials Vn and Vp in Eq. ( 14 ) turn out to be rapidly growing functions 



of density, and one can anticipate that eventually reactions such as 

p + e^ A + Ve and/or n + ^ + 17^ (16) 

may produce hyperons. Hyperons can appear once the corresponding /3-equi- 
librium conditions are satisfied, i.e., /i„ — or/and fj,n + Pe = Py,- ■ At the 
threshold, where pfk = or ppy,- — 0, one can expect that |Va| ^ toa and 
|Vx;-| ^ m-s- and thus /^a — 'tia and pL^- ~ m^- . Since mA and mj^- are 
larger than the nucleon mass by only about 200 MeV these hyperon^ are good 
candidates for an "exotic" form of matter in neutron stars. Along similar lines, 
the lightest mesons, pions and/or kaons, may also appear stably, and can form 
meson condensates. At even larger densities, the ground state of matter is 
likely to be one of deconfined quarks. All these possibilities depend crucially 
on the strong interactions terms, Vn and Vp. Figure [3] illustrates our present 
understanding (or misunderstanding) of the interior of a neutron star, with a 
black question mark "?" in its densest part. The outer part of the star, its crust, 
is described briefly in the following subsection. 

When only nucleons, plus leptons as implied by charge neutrality and con- 
strained by /3-equilibrium, are considered, the EOS can be calculated with much 
more confidence than in the presence of "exotic" forms of matter. For illustra- 
tive puposes, we will generally employ the EOS of Akmal, Pandharipande & 
Ravenhall [16] ("APR" hereafter) in presenting our results. 

Relativistic expressions for fin and /jp also exist, but are omitted here in the interest of 
simplicity. 

*The S"*" is less favored as its /3-equilibrium condition is = /^p = — Me. Heavier 

baryon are even less favored, but cannot a priori be excluded. 
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Although there is no evidence that any observed neutron star or pulsar might 
actually instead be a pure quark star, theory allows this possibility. Such a star 
would be nearly completely composed of a mixture of up, down and strange 
quarks, and would differ from a neutron star in that it would be self-bound 
rather than held together by gravity. 

The reader can find a more detailed presentation and entries to the key 
literature in [15j . 

3.2 The Neutron Star Envelope and Crust 

In the outer part of the star, where p ^ Pnuc/2, a homogeneous liquid of nucle- 
ons is mechanically unstable (known as the spinodal instability). Stability is, 
however, restored by the formation of nuclei, or nuclear clusters. This region, 
called the crust, has a thickness of 1 km. 

Above the surface, where the pressure approximately vanishes, we might 
expect the presence of an atmosphere, but there is the possibility of having 
a solid surface, condensed by a sufficiently-strong magnetic field [17]. A few 
meters below the surface, ions are completely pressure-ionized (the radius of 
the first Bohr orbital is larger than the inter-nuclear distance when p ^ 10'* g 
cm~^). Matter then consists of a gas/liquid of nuclei immersed in a quantum 
liquid of electrons. When p w 10^ g cm^'^, fie is of the order of 1 MeV and the 
electrons become relativistic. Here, and at higher densities. Coulomb corrections 
are negligible - electrons form an almost perfect Fermi gas. However, Coulomb 
corrections to the ions are not negligible. From a gaseous state at the surface, 
ions will progressively go through a liquid state (sometimes called the ocean) 
and finally crystallize, at densities between 10^ up to ^ 10^" g cm~^ depending 
on the temperature (within the range of temperatures for which neutron stars 
are thermally detectable). These low-density layers are commonly referred to 
as the envelope. 

With growing p, and the accompanying growth of /i^, it is energetically 
favorable to absorb electrons into neutrons and, hence, nuclei become progres- 
sively neutron-rich. When p ^ Ax 10*^ g cm""^ (the exact value depends on the 
assumed chemical composition), one achieves the neutron drip point at which 
the neutron density is so much larger than that of the protons that some neu- 
trons become unbound (i.e., pn > 0). Matter then consists of a crystal of 
nuclei immersed in a Fermi gas of electrons and a quantum liquid of dripped 
neutrons. The region containing dripped neutrons is usually called the inner 
crust. In most of this inner crust, because of the long-range attractive nature 
of the nucleon-nucleon interaction, the dripped neutrons are predicted to form 
a superfluid (in a spin-singlet, zero orbital angular momentum, state ^Sq). All 
neutron stars we observe as pulsars are rotating. While a superfluid cannot 
undergo rigid body rotation, it can simulate it by forming an array of vortices 
(in the cores of which superfluidity is destroyed). (See, e.g., [18].) The resulting 
structure is illustrated in inset A of Fig. |3] 

At not too-high densities, nucleons are correlated at short distances by strong 
interaction and anti-correlated at larger distances by Coulomb repulsion be- 
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tween the nuclei, the former producing spherical nuclei and the latter resulting 
in the crystallization of the matter. As p approaches '-^ O.OSpnuc, the shapes 
of nuclei can undergo drastic changes: the nuclear attraction and Coulomb re- 
pulsion length-scales become comparable and the system is "frustrated" . From 
spherical shapes, as the density is increased, nuclei are expected to deform, be- 
come elongated into 2D structures ("spaghetti"), and then form ID structures 
("lasagna"), always with denser nuclear matter surrounded by the dilute neu- 
tron gas/superfluid which occupies an increasing portion of the volume. When 
the phases achieve approximately equal volume fractions, the geometry can in- 
vert, with dripped neutrons confined into 2D ("anti-spaghetti" or "ziti") and 
finally 3D ("swiss cheese") bubbles. The homogeneous phase, i.e. the core of 
the star, is reached when p ~ 0.5 — O.Gpnuc- This "pasta" regime is illustrated in 
inset B of Fig. [s] and is thought to resemble a liquid crystal [19] . A compilation 
of the most recent progress on neutron star crust physics can be found in the 
book dD]. 



4 Pairing in the Neutron Star Interior 

4.1 General Considerations 
Expectations from measured phase shifts 

As a two-particle bound state, the Cooper pair can appear in many spin-orbital 
angular momentum states (see the left panel of Fig. |4|. In terrestrial super- 
conducting metals, the Cooper pairs are generally in the ^Sq channel, i.e., spin- 
singlets with L = orbital angular momentum, whereas in liquid ^He they 
are in spin-triplet states. What can we expect in a neutron star? In the right 
panel of Fig. [4] we adapt a figure from one of the first works to study neutron 
pairing in the neutron star core |21j showing laboratory-measured phase-shifts 
from nucleon-nucleon scattering. A positive phase-shift implies an attractive 
interaction. From this figure, one can expect that nucleons could pair in a spin- 
singlet state, ^Sq, at low densities, whereas a spin-triplet, '^P2, pairing should 
occur at high densities. We emphasise that this is only a presumption (phase 
shifts reflect free-space interaction) as medium effects can strongly affect particle 
interactions. 



The energy gap 

In a normal Fermi system at T = 0, all particles are in states with energies 
e < ep. When T > 0, states with energies e ^ cp can be occupied (left panel 
of Fig. [5]) resulting in a smearing of the particle distribution around in a 
range ^ ksT . It is precisely this smooth smearing of energies around which 



produces the linear T dependence of c„. Sec. 7.1 and the or dependence 
of the neutrino emissivities, Sec. |6.5[ 

In a superfluid/superconducting Fermi system at T = 0, all particles are in 
states with energies e < ep (actually, e < ep — A). For nonzero temperatures 
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Figure 4: Left panel: Possible spin-angular momentum combinations for 
Cooper-pairs. Right panel: Phase shifts for N-N scattering as a function of 
the laboratory energy (middle axis) or the neutron Fermi energy and density 
for a neutron star interior (lower axis). Adapted from |21j . 



that permit the presence of Cooper pairs (and hence a gap A(r)), states with 
energy e > e^? + A can be populated. However, in contrast to the smooth filling 
of levels above e^? in the case of a normal Fermi liquid, the presence of the 
2A(r) gap in the spectrum implies that the occupation probability is strongly 
suppressed by a Boltzmann-likc factor ^ cxp[— 2A(T)/fcs2^]- As a result, both 
the specific heat of paired particles and the neutrino emissivity of all processes 
in which they participate are strongly reduced. 

The phase transition 

The transition to the superfluid/superconducting state through pairing d la 
BCS is usually a second order phase transition and the gap A(r) is its order 
parameter (see central panel of Fig. [6]). Explicitly, A(T) = when T > Tc, 
the critical temperature, and, when T drops below T^, A(T) grows rapidly but 
continuously, with a discontinuity in its slope a,t T — Tc- There is no latent heat 
but a discontinuity in specific heat. (Examples: superfluid O normal fluid; fer- 
romagnetic O paramagnetic.) In the BCS theory, which remains approximately 
valid for nucleons, the relationship between the zero temperature gap and Tc is 

A(T = 0) ~ 1.75 fcsT,. (17) 

In a first order phase transition there is a discontinuous change of A(T) at Tc 
and the transition occurs entirely at Tc (see left panel of Fig.[6|. There is a latent 
heat due to the entropy difference between the two states. (Examples: solid O 
liquid; liquid ^ gas below the critical point.) In a smooth state transition there 
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Nomal Fermi Liquid Superfluid Fermions 



e 




Figure 5: Comparison of quasi-particle spectra, e(fc), for a normal and a super- 
fluid Fermi liquid. The reorganization of particles at e ^ ep into Cooper pairs 
results in the development a gap 2A in the spectrum so that no particle can 
have an energy between ~ A and + A. 



First Order Phase Transition 




Second Order Phase Transition 
A" 



Smooth State Transition 





Figure 6: Temperature evolution of the state of a system parametrized by an 
"order" parameter, A(T). 



is a continuous change of A(r) with no critical temperature (see the right panel 
of Fig. [6|. (Examples: liquid -;-> gas above the critical point; atomic gas o 
plasma.) 

A simple example 

A simple model can illustrate the difhculty in calculating pairing gaps. Con- 
sider a dilute Fermi gas with a weak, attractive, interaction potential U. The 
interaction is then entirely described by the corresponding scattering lengtlj^, a, 
which is negative for an attractive potential. In this case, the model has a single 
dimensionless parameter, |a|/cF, and the dilute gas corresponds to \a\kF ^ 1- 
Assuming the pairing interaction is just the bare interaction U (which is called 
the BCS approximation), the gap equation at T = can be solved analytically, 

^The scattering length a Is related to U hy a = {in / inh^ JUo with Uy^ = J cPr exp(ik • 
r) C/(r). 



14 



giving the weak- coupling BCS- approximation gap: 

L \ r TT 

(18) 



A(fcp) n , ^Bcsikp) = ^ ^ ) exp 



2\a\ki 



This result is bad news: the gap depends exponentially on the pairing potential 
U. The Cooper pairs have a size of the order of ^ ~ Tivp/^ (the coherence 
length) and thus S^kp oc exp[+7r/2|a|/ci?] 3> 1 in the weak coupling limit. There 
is, hence, an exponentially growing number of other particles within the pair's 
coherence length when |a|fci^- — > 0. These particles will react and can screen, or 
un-screen, the interaction. Including this medium polarization on the pairing is 
called beyond BCS, and in the weak coupling limit its effect can be calculated 
analytically j22j, giving 

A(fc^) AGMBikp) = j^ABcsikp) - OAdABcsikp) (19) 

So, screening by the medium can reduce the gap by more than a factor two, 
even in an extremely dilute system. 

4.2 Calculations of Pairing Gaps 

A significant amount of work has been devoted to the calculation of pairing 
gaps in the neutron star environment: see, e.g., [531 IMl Hi] or A. Schwenk's 
contribution to this volume for reviews. Below we first briefly describe the 
Gorkov formalism [261 EH El] that will allow us to set up the stage for the 
presentation of representative results for nuclcon gaps. Specifically, we will 
address the neutron -"^Sq and '^P-F2 and the proton ^Sq and ■^P-F2 gaps and 
briefly mention hyperon gaps. The effects of pairing on the thermal evolution 



of neutron stars are described in Sec. 4.3 Calculations of pairing gaps in quark 
matter will be described in Sec. [U 

General formalism 

Several significant effects of pairing are due to the change in the quasi-particle 
spectrum that is obtained from the poles of the propagator G. In Fig.[7j Eq. (A), 
we show the definition of G, the two point Green's function with one particle in 
and one particle out, and its Dyson equation which relates the free propagator 
and the self energy E in the case of a normal system. The resulting quasi-particle 
spectrum is then 

u =e{k, Lj) = —— + S(fc, u)-£F^ -^{k - kp) . (20) 
2m TO* 

In obtaining the right-most result above, we assumed the system to be isotropic 
and the spectrum is evaluated for k ^ kp with the effective mass m* defined 
through 

de{k,uj = 0) 



dk 



fi^kp 



(21) 



15 



(A) zG(l,2) = (T{^(l)V't(2)}) = = ^ + ^-©-^ 

Figure 7: The normal state propagator and its Dyson equation. 



The resulting spectrum ep + ^jJ is depicted in the left panel of Fig. [5) 

In the presence of a pairing instability, and the concomitant development of 
a condensate, an anomalous propagator F and its adjoint can be defined, 
see Eq. (B) and Eq. (B') in Fig. [sj with their corresponding Gorkov equations 
that replace the Dyson equation. In addition to the self energy E, the Gorkov 
equations feature an anomalous self energy, or the gap function A. The prop- 
agator F violates particle number conservation as it propagates a hole into a 
particle, and vanishes in the absence of a condensate in a normal system. The 
gap function is a 2 x 2 matrix in spin space 

A(fc,-)=f ) • (22) 

In the case the ground state is assumed to be time-reversal invariant, A(fc,a;) 
(denoted A for short below) has a unitary structure satisfying 

AA^ = A^A = A^ i (23) 

where 1 is a 2 x 2 unit matrix, and 

A^ = A2(fc,a;) =detA(fe,w). (24) 

The quantity A above will appear as the energy gap in the quasi-particle spec- 
trum. 

The normal propagator G is also modified, as depicted in Eq. (C) in Fig. [sj 
Solving the Gorkov equations gives G — {uj + e) / D with 

D(fc,a;) =a;2-e(fc,a;)2- A(fc,w)2 (25) 
and its modified poles yield a quasi-particle spectrum with two branches: 



w = ±Ve(fe,w)2 + A(fe,w)2. (26) 

The resulting spectrum e^+i^ is depicted in the right panel of Fig. [5] Note that, 
G and S are also 2x2 matrices in spin space both being diagonal in structure. 
The solution of the Gorkov equation for F gives F ^ KyO. Finally, E and A 
are defined by Eq. (D) and (E) in Fig. [sjfrom a kernel K. Equation (D) is the 
gap equation and reads 

A ■ [ ^ (fca,-fc/3|X|fcV,-fc'^')A^,^,(fcO 

where k = (k,uj), k' = {k',uj') and a, /3, ... denote spin indices. 
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(B) ^F(1,2) = 








(B') iFHl,2) 


= (r{V't(i)^t(2)}) = 


— + - 




(C) zG(l,2) = 











(D)H|)-=^J^ (E)^^= IK 

Figure 8: The Gorkov equations 

Solving the full set of equations in Fig. [8] requires many approximations, no- 
tably in the choice of the kernel K. In particular, different approximate kernels 
K are used in (D) and (E) as it appears in the particle- hole channel in (D), 
whereas in (E) it is in the particle-particle channel. In the BCS approxima- 
tion, the self energy is calculated in the normal phase (e.g., with a Briickner- 
Hartree-Fock or BHF scheme) with its energy dependence being neglected, i.e. 
'E{k,uj) — )■ S(fc), and the kernel for A is simply taken as the bare interaction. 
The latter can be a two-body force (2BF) from a potential adjusted to labora- 
tory N-N scattering data or that derived from an effective interaction such as 
a Skyrme force. Recently, the effect of three-body forces (3BF), absent in the 
laboratory N-N scattering experiment, has been considered. Inclusion of 3BF 
is necessary to reproduce the nuclear saturation density; they are, in the bulk, 
repulsive and their importance grows with increasing density. Even at the 2BF 
level, a severe problem is encountered: to date, none of the N-N interaction 
models reproduce the measured phase shifts for E\ah larger than 300 MeV in 
the channels needed (particularly for ^P2) for the conditions prevailing in the 
core of a neutron star. The source of this problem is easy to identify; beyond 
290 MeV (the threshold for pion production) , inelastic channels begin to become 
dominant. 

Models beyond BCS have proceeded in two directions. In the first approach, 
in the gap equation, the bare interaction is supplemented by the inclusion of 
short-range correlations. In a further step, long-range correlations to account for 
medium polarization are added. As illustrated above with the weak-coupling re- 



sult of Eq. ( 19 1, polarizations effects can be significant. In the second approach, 
which goes beyond the BHF level, the self energy E is calculated by including 
its energy dependence. Calculations of E in the paired phase, see Eq. (D) in 
Fig. [Sj do not yet exist and are necessary. In the following, we will neglect the 
energy dependence of the gap, i.e., write it as A(fc) instead of A{k,uj). 

Pairing in single spin-angular momentum channels 

In pairing calculations, the potential and the gap function are usually expanded 
in partial waves so as to focus on specific spin-angular momentum channels, A = 
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(s, j). At low kpn, 01 kpp, it is theoretically predicted that the preferred channel 
is A = (0,0) in ^-wave, i.e., the spin-singlet ^Sq. At large Fermi momenta, the 
^Sq interaction becomes repulsive and the preferred channel is A = (1,2) in P 
and F waves (the mixing being due to the tensor interaction [28j ) , i.e., the spin- 
triplet ^P-F2. In the ^Sq channel, which has also been called the "A" phase, the 
gap is spherically symmetric and can be written as 

A phase (iSo) : A(o.o)(fe) = (^_^^^^ ^^^^^ [Aikp) = energy gap] (28) 

In the ■^P-F2 channel, A;^ has contributions from all possible orbital angular 
momenta I and their rrij components, i.e., Aa = J2i m ^Tx i^)^Tx i^)' where 
the (^) ^'"'3 2x2 spin matrices describing the angular dependence of A;^ 
which is thus not spherically symmetric. Microscopic calculations restricted to 
the '^P2 channel [221 130] indicate that the largest component of Ax corresponds 
to the rrij = sub-channel or, possibly, the rrij = ±2 channels, sometimes called 
the "B" and "C" phases, respectively. For these two special cases, the energy 
gap A{kF) is given by [29] 

B phase (3p2,m, = 0) : A'{kp) = 1 [A^^xikpf ^^^g^^^' ^^^^ 

Cphase(3p2,m, =±2): A'{kp) = [Alxikp)]^ ^^"^ . (30) 

where 9k is the angle between kp and the arbitrary quantization axis. Notice 
that in the B phase, the gap is nodeless whereas in the C phase it vanishes on 
the equator, = tt/2. 

Temperature dependence of A and 

The preceding discussion was restricted to the zero temperature case. It is natu- 
rally extended to finite temperature whence the gap becomes A(fc; T). However, 
effects of thermal excitations are important only for values of fc ~ fc^- We will 
often omit either of the arguments kp or T when not necessary, but they are 
always implied (as is its uj dependence). Notice that microscopic calculations 
are often limited to the T = case only. 

The relationship between the critical temperature Tc for the phase transition 
and the energy gap A{kp) is approximately given by the usual result 

fcsTcW 0.57 A(fci.;r = 0) (31) 

for all three phases A, B, and C [211 [30], where A{kp;T) is obtained by angle 
averaging of A'^{kp,T) over the Fermi surface 

[Aikp;Tf^ JJ^A'ikp;T). (32) 

Obviously, A{kp;T) — A{kp;T) for an isotropic ^Sq gap. The temperature 
dependence of the energy gap A{kp,T) for ^Sq pairing and of the angle averaged 
A{kp; T) for the '^P2 pairing in the rrij — case are shown in Fig. fel 
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T/T, 



Figure 9: Temperature dependence of the energy gap A{kp, T) for ^Sq pairing 
and of the angle averaged A{kp;T) for the ■^P2 pairing in the mj — case. 
Continuous hues show the commonly used analytical fits of [31 whereas the 
dots are from the calculations of [5D1. '^P2 pairing in the rrij — 2 case results in 
values very close to the mj = case. 



The isotropic ^Sq neutron gap 



In Fig. 10 we show sets of predicted Tc for the neutron ^Sq pairing in uniform 



pure neutron matter. The two dotted lines marked "BCS" and "GMB" show 



the simple analytical results of Eq. (18 1 and Eq. (19), respectively, with a = 



— 18.5 fm and = 0.56 A{kp). Formally, these results are only valid when 
\a\kF < 1, i.e., kp < 0.1 frn'^ The curve "SCLBL", from [35], illustrates 
the results of a numerical solution of the gap equation, using the Argone Vn 
N-N potential, within the BCS approximation. The results merge with the 
"BCS" curve in the weak coupling limit kp —i' 0. Also shown are results from 
calculations that take into account more sophisticated medium effects, including 
medium polarization (with different schemes): "CCDK", from [33], employed 
a variational method within the correlated basis functions scheme, "WAP", 
from |34j . employed an extension of the induced interaction scheme, whereas 
"SFB", from [35], went beyond [34] with renormalization group methods. In 



line with the simple GMB result of Eq. ( 19 1, these model calculations show that 



polarization has a screening effect that quenches the gap, by a factor ~ 3. These 
three calculations yield some agreement, particularly in the predicted maximum 
value of Tc, but with a non- negligible difference in the density dependence. The 
other two curves show more recent results: "GIPSF", from [36], utilizes the 
auxiliary-field diffusion Monte Carlo technique while "GC", from [37], stems 
from a Quantum Monte Carlo calculation. These last two models result in gaps 
that are intermediate between the previous models and the BCS approximation; 
moreover, they converge toward the GMB value when kp ^ 0. 
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Figure 10: Some theoretical predictions of and A, vs neutron kp, for the 
neutron ^Sq gap in uniform pure neutron matter. The value of kp corresponding 
to the transition from the crust to the core is indicated. See text for description. 



Transcribed to the neutron star context, the range of Fermi momenta for 
which these neutron -'^Sq gaps are non vanishing corresponds mostly to the 
dripped neutrons in the inner crust. The presence of nuclei, or nuclear clus- 
ters in the pasta phase, may modify the sizes of these gaps from their values in 
uniform matter. The coherence length ^ of the dripped neutrons is larger than 
the sizes of nuclei, leading to proximity effects. This issue has received some 
attention, see, e.g., [35113211101111], and position dependent gaps, from inside 
to outside of nuclei, have been calculated. However, in most of the crust ^ is 
smaller than the internuclear distance, and the size of the gap far outside the 
nuclei is close to its value in uniform matter. 



The isotropic ^Sq proton gap 

The magnitudes of proton ^Sq gaps are similar to those of neutrons, but with the 
important difference that, in the neutron star context in which beta equilibrium 
prevails, protons are immersed within the neutron liquid, and constitute only 
a small fraction of the total baryon number (3 to 20% in the density range 
where they are expected to be superconducting). Proton-neutron correlations 
cause the effective mass of the proton to be smaller than that of the neutron, 
a simple effect that reduces the size of the proton ^Sq gap compared to that of 
the neutron. 

Several theoretical predictions of Tc for the proton ^Sg gap are shown in 



Fig.[n] "CCY" from [H], "T" from [13], and "AO" from [H] that are among 
the first historical calculations, whereas "BCLL" from [30], "CCDK" from [33] . 
and "EEHO" from [IS] are more recent results. All of these calculations were 
performed within the BCS approximation and very few works have gone beyond 
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Figure 11: Some theoretical predictions of and A, vs proton kp, for the 
proton ^Sq gap in /3-equiHbrium uniform neutron-proton matter. The value of 
kp corresponding to the transition from the crust to the core is indicated: values 
on the right of this line correspond to the neutron star core but values on the 
left are not realized since protons in the crust are confined within nuclei which 
are finite size systems while this figure presents results for infinite matter. On 
the top margin are marked the values of the proton kp at the center of a 1.0, 
1.4, 1.8, and 2.0 Mq star built with the APR EOS [H]. See text for description. 

BCS for the proton ^Sq gap. Among the latter, we show results from jl^ : 
these authors used either only two body forces in the interaction kernel, curve 
"BS2bf" , or two body forces supplemented by the inclusion of three body forces, 
curve "BS2BF+3BF" which shows that three body forces are repulsive in the ^Sq 
channel. These "BS" results also include effects of medium polarization. Recall 
that for the ^Sq pairing of neutrons in pure neutron matter, polarization has a 
screening effect and quenches the gap. However, in neutron star matter, where 
the medium consists mostly of neutrons, the strong np-correlations result in 
medium polarization inducing anti- screening [47j for the ^Sq pairing of protons. 

The anisotropic '^P-F2 neutron (and proton) gap 

The ^So neutron gap vanishes at densities close to the crust-core transition 
and the dominant pairing for neutrons in the core occurs in the mixed '^P-F2 
channel. Uncertainties in the actual size and the range of density in which this 
gap persists are however considerable. As previously mentioned, a major source 
of uncertainty is the fact that even the best models of the N-N interaction in 
vacuum fail to reproduce the measured phase shift in the ■^P2 channel [48]. Also 
significant are the effects of the medium on the kernel and 3BF, even at the 
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level of the BCS approximation. It was found in [35] that 3BF at the Fermi 
surface are strongly attractive in the '^P-F2 channel in spite of being repulsive in 
the bulk. Moreover, due to medium polarization a long-wavelength tensor force 
appears that is not present in the interaction in vacuum and results in a strong 
suppression of the gap [SU] . 
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Figure 12: Left panel: some theoretical predictions of for the neutron '^P-F2 
gap in uniform pure neutron and /3-equilibrium matter. See text for description. 
Right panel: some phenomenological models of Tc for the neutron ■^P-F2 gap 
used in neutron star cooling simulations. Models "a", "b", and "c" are from 
[5T] and [52], model "a2" from ,53J. On the top margin are marked the values 
of kpn at the center of a 1.0, 1.4, 1.8, and 2.0 Mq star buih with the APR EOS 

m- 



Figure 12 shows examples of theoretical predictions of Tc for the neutron 
^P-F2 gap. The three dotted lines show some of the first published models: 
"HGRR" from [53], "T" from and "AO" from [29 . The four continuous 
lines show results of models from [48, calculated using the Nijmegen II ("Nijll"), 
Nijmegen I ("NijI"), CD-Bonn ("CDB"), and Argonne Vis ("AV18") potentials 
(displayed values are taken from the middle panel of Figure 4 of pSj). The 
results of these four models start to diverge at kpn above 1.8 fm^^ and illustrate 
the failure of all four N-N interactions models to fit the ^P2 laboratory phase- 
shifts above E'lab — 300 MeV. All of these calculations were performed for pure 
neutron matter using the BCS approximation. 

In the case of the "'^Sq gap, medium polarization is known to result in screen- 
ing and to reduce the size of the gap. In the case of a ^P2 gap, polarization 
with central forces is expected to result in anti-screening and to increase the 
size of the gap. However, Schwenk & Friman |50j showed that spin-dependent 
non-central forces do the opposite and strongly screen the coupling in the '^P2 
channel, resulting in a Tc lower than 10^ K: this "SF" value is indicated in the 
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Figure by an arrow. 

The two dashed hnes in Fig. 12 present resuhs from where the "2BF" 
model only considers 2-body forces (from the Argonne Vis) while the "3BF" 
model includes a meson exchange model 3-body force: the result is a growing 
^P-F2 gap which shows no tendency to saturate at high density. This work, 
for /3-equilibrium matter, moreover emphasized the importance of the proton 
component. 

Other delicate issues are the effect of the proton contaminant and the likely 
development of a tt" condensat^ which also strongly affects the size of the 
neutron (and proton) gap(s). 

In summary, the size and extent in density of the neutron ■^P-F2 gap in the 
neutron star core is poorly known. Given these large uncertainties in the size 
of the neutron '^P-F2 gap (about three orders of magnitude) and the fact that 
neutrino emissivity is suppressed by an exponential Boltzmann-like factor, this 
gap is often considered as a free parameter in neutron star cooling models. The 
extreme sensitivity of the cooling history on the size of this gap can be utilized 
to one's advantage by inverting the problem, as it may allows us to measure it by 



fitting theoretical models to observational data |57j. The right panel of Fig. 12 
presents the phenomenological neutron '^P-F2 gaps used in cooling calculations 
in a later section. 

In the case of protons, their '^P-F2 gaps have generally been overlooked due to 
their small effective masses, and considered to be likely negligible [5H]. However, 
in view of the strong enhancing effect of the 3BF on the neutron ^P-F2 gap, this 
issue has been reconsidered in [49] where it was shown that the proton ■^P-F2 
gap can be sizable. 



Hyperon gaps and nucleon gaps in hyperonic matter 

Many calculations of dense matter indicate that strangeness-bearing hyperons 
will be present in neutron star matter once the neutron chemical potential ex- 



ceeds the rest masses of hyperons [59] (see Sec. 3.1). In the likely presence of 
hyperons (denoted by Y) arise the issues of, first, the effect of their presence on 
the nucleon gaps, and second, the possibility of hyperon pairing. Nucleon gaps 
in the presence of hyperons have been studied in [Aff and [5(T: depending on 
the N-N, N-Y, and Y-Y interaction models employed, the nucleon gaps may be 
either enhanced or reduced by the presence of hyperons. 



Since the suggestion in and the first detailed work of [55], hyperon gaps 
have attracted some attention. All uncertainties present in the nucleon case 
immediately translate to the hyperon case. An additional problem is that very 
little is known about hyperon-hyperon interaction^ which is generally guessed 

^In the presence of a charged n~ condensate a new Urea neutrino emission pathway is 
open, see Table ^ The development of a neutral tt" condensate has, however, little effect on 
neutrino emission. 

'^Some experimental information is available from hypernuclei |63l I64| and hadronic 
atoms | 65| . but the data do not yet uniquely determine the hyperon- nucleon or hyperon- 
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from theory by extrapolation from the N-N interaction. At densities not much 
larger than their threshold densities for appearance, hyperons have low enough 
concentration that they will pair in the ^Sq channel. In both cases of A's and 
S~'s, the estimated gap sizes are similar to those of nucleons. We refer the 
reader to [171 |SH1 El] and references therein for details. 

At very high densities, the presence of deconfined quark matter is also likely. 
Quarks are expected to pair and form a color superconductor. This subject has 
developed into a field of its own, and we dedicate Sec. [5] for a brief account and 
refer the reader to other articles in this monograph. 



4.3 Effects of Pairing for Neutron Star Cooling 

The occurrence of pairing leads to three important effects of relevance to neutron 
star cooling: 

A) Alteration, and possible strong suppression when T <^ Tc, oi the specific 

heat Cv of the paired component. 

B) Reduction, and possible strong suppression when T <^ Tc, of the emissivity 

of the neutrino processes the paired component is involved in. 

C) Triggering of the "Cooper pair breaking and formation" (PBF) with con- 

comitant neutrino pair emission which is very efficient in the case of spin- 
triplet pairing. 

These effects are direct consequences of the development of the energy gap A(fc) 



and the resulting two branches in the quasi-particle spectrum, Eq. ( 26 1 . The 
gap severely limits the available phase space when T <^ Tc and the spectrum is 
usually treated in the effective mass approximation with an angle averaged gap. 



see Eq. (32) 



e{k) = ±^ [hvF{k - kpW + A(fcj., T) , (33) 

where vp = hkp/m* . 

In cooling calculations, these effects are introduced through "control func- 
tions" : 

Cv Rc cv and Rx . (34) 

There is a large family of such control functions for the various types of pair- 
ing and the numerous neutrino processes "X" . For nodeless gaps, the i?'s are 
Boltzmann-like factors ~ exp[—2D(T)/kBT] and result in a strong suppression 
when T <ti Tc, whereas for gaps with nodes the suppression is much milder. 
Regarding the specific heat, there is a sudden increase, by a factor ~ 2.4 at 
T = Tc, followed by a reduction at lower T. Examples of such control functions 
are shown in Fig. [Ts) 



hyperon interaction. Future work in lattice QCD |66| may prove fruitful, but current results 
are limited to unpliysically large pion masses. 
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Figure 13: Control functions for neutrino emission from the modified Urea pro- 
cess (as, e.g., n + n n + n + vV) (left panel) and the specific hea t (r ight panel), 
in the presence of ^Sq pairing and ■^P2 in the B phase (see Eq. (29)), from the 
analytical fits of [31] and [70] . 



The effect C, neutrino emission from the formation and breaking of Cooper 
pairs [711 172] . can be interpreted as an inter-band transition (as, e.g., n — > 
n + vV) where a neutron/proton quasiparticle from the upper (-I-) branch of 
the spectrum of Eq. (33) falls into a hole in the lower (— ) branch. Such a 
reaction is kinematically forbidden by the excitation spectrum of the normal 
phase, Eq. (20), but becomes possible in the presence of an energy-gap, Eq. (33). 



This process is described in more detail in Sec. |6.6[ The resulting emissivity 
can be significantly larger than that of the modified Urea process (as, e.g., 
n + rt— >n + rt + vV) in the case of spin-triplet pairing. 



5 Superfluidity in Dense Quark Matter 

The central densities of neutron stars can exceed the nuclear density Pnuc ~ 
2.7 X 10^^ g cm~'^ by significant amounts. At sufficiently high densities, a de- 
scription of neutron star interiors in terms of nucleons becomes untenable and 
sub-nucleonic degrees of freedom, namely quarks, must be invoked. Interac- 
tions between quarks is fundamentally grounded in Quantum Chromodynamics 
(QCD), the theory of strong nuclear interactions. The theory has a gauge sym- 
metry based on the Lie group SU(3), and the associated charge is referred to 
as "color". QCD is asymptotically free: interactions between quarks mediated 
by gluons become weak at short distances, or equivalently, high densities. At 
low densities, strong interactions "confine" quarks into neutrons and protons 
which are color neutral. Asymptotic freedom guarantees that, at some large 
density, the ground state of zero-temperature matter will consist of nearly-free, 
"deconfined" quarks [75] . 

QCD has been amply tested by experiments at high energies where asymp- 
totic freedom has been confirmed 74]. Lattice-gauge calculations of hadron 
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masses, and of baryon-free matter at finite temperature, liave made enormous 
strides in recent years [7S]. However, first-principle QCD calculations for fi- 
nite baryon density have been stymied due to the fermion-sign problem in lat- 
tice gauge calcuations. While many guesses are available, perturbation theory 
is unable to accurately predict the density at which the deconfinement phase 
transition occurs. However, it is possible that the phase transition occurs at 
a density lower than the central density of some (or even all) neutron stars. 
In that case, at least some neutron stars will contain deconfined quark matter. 
Such objects are referred to as hybrid quark stars. In their cores, up (u), down 
(d) , and strange (s) quarks are the principal degrees of freedom, the other three 
quark flavors (charm (c), bottom (6), and top (t)) being excluded because of 
their large masses. It is theoretically possible that the energy of zero-pressure 
strange quark matter has a lower energy than Fe [76l [77l [78j , in which case a 
hybrid star would be metastable or unstable, and nucleonic matter would spon- 
taneously convert into strange quark matter, creating pure quark stars that are 
self bound (TnHHU]. There is no experimental or observational evidence for pure 
quark stars, however, and we do not consider them further in this contribution. 

5.1 Pairing in Quark Matter 

Cooper pairing between quarks was first investigated in the late 1970s [STl 15^ . 
As all quarks are charged, pairing between quarks is often referred to as color 
superconductivity because the paired phase breaks the SU(3) gauge symmetry 
of QCD. First estimates of pairing gaps in quark matter were of order 1 MeV. In 
this case, neutron stars with deconfined quark matter cool very rapidly through 
the quark direct Urea processes (see Table [T]). Because some older neutron stars 
are observed to be relatively warm, this naturally implies that not all neutron 
stars can contain quark matter. 

This situation changed drastically with the discovery ^S31 HI] that color 
superconductivity implies gaps as large as 100 MeV (see a recent review in 
[55]). These works suggested two possibilities: either the so-called "color-flavor- 
locked" (CFL) phase in which all nine combinations of flavor (up, down, strange) 
and color (red, green, blue) participate in pairing, or the "2SC" phase where 
only four of the nine combinations pair (corresponding to up and down quarks 
which are either red or green). 

5.2 Theoretical Descriptions of Dense Quark Matter 

There are several formalisms which have been applied to describe color super- 
conducting quark matter. High-density effective theories (HDET) were first 
developed in 1990s j86j [87] and then developed further for color superconduc- 
tivity in [SHI |89] . To construct an effective theory, one begins with the QCD 
Lagrangian, rewrites it in terms of a 1//^ expansion (/x being the chemical poten- 
tial), and then integrates out hard gluons and fermionic modes corresponding 
to the Dirac sea. 
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Another commonly used alternative consists of using Nambu-Jona-Lasinio 
(NJL) 9QJ models. The original NJL model was a theory of strong interactions 
before the advent of QCD. The four-fermion interaction of the NJL Lagrangian 
bears close resemblance to that in the BCS theory of superconductivity and 
gives rise to analogous effects. Originally framed in terms of nucleon fields, 
quartic interactions serve to give the nucleon its mass through a self-energy 
generated by the formation of a condensate. Modern versions involve quark 
fields that develop a "quark condensate" which is then related to the mass of 
constituent quarks [91]. NJL models for color superconductivity presume that 
gluonic degrees of freedom have been integrated out resulting in point-like cou- 
plings between quarks. For a review of the NJL model applied to dense quark 
matter see |92| . Confinement is sometimes implemented by the addition of 
Polyakov loop terms giving rise to "PNJL" models. HDET and NJL methods 
give qualitatively similar results, but the NJL Lagrangian is a bit more trans- 
parent, so we describe some of its details here. A chiral SU(3) Lagrangian with 
superconducting quarks (adapted from [93] ) is 

8 

+G5^[(9A}(7)' + (g»75A;g)' 

a=0 

+ Ga^^ [qtai-iiiki-afiiq^f^ {^'a''^i'j'k'f-a'l3'-i'qrf}') 
k 7 

+ G/^^^{qtailb^ijk^al3iq%) {qt'a'il5^i'3'k'^a'l3'i'qrfi') , (35) 
k 7 

where Roman indices are for flavor and greek indices are for color, except for 
a which enumerates the SU(3) matrices, rriij is the quark mass matrix, fJ-ij.ap 
is the chemical potential matrix, q is the quark field, = Clf" , and e is the 
Levi-Civita tensor. The first term is the Dirac Lagrangian which describes free 
relativistic massive quarks at finite density. The second term is a combination of 
four quark fields which model non-superfluid quark-quark interactions and obeys 
the SU{3)l X SU{3)u chiral symmetry present in QCD. The third and fourth 
terms, which give rise to color superconductivity, are the chirally symmetric 
analog of the second term in the quark-quark channel. QCD breaks C/a(1) 
symmetry, and so these four-fermion interactions can be supplanted by six- 
fermion interactions in order to do the same [M] . 

The first step in obtaining the thermodynamic potential in the mean-field 
approximation is to replace the quark bilinears qiQi and qiai'^h^'^^'' ^"^^^ qfp with 
their ground-state expectation values. The former is the quark condensate as- 
sociated with the breaking of chiral symmetry and the latter gives rise to the 
superconducting gap, A'^'*'. Having made this replacement, the non-constant 
terms in the Lagrangian take the form qMq^ where M is a matrix represent- 
ing the inverse propagator. This matrix can be diagonalized in the standard 
way to obtain the individual quark dispersion relations and the thermodynamic 
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(36) 



where give the energy eigenvalues and i runs over 3 flavors, 3 colors, and 
the Dirac indices (36 total). In this formulation, there is a manifest parallelism 
between the quark condensates, (gg), and the superconducting gaps. The mini- 
mum of the thermodynamic potential with respect to the superconducting gap 
gives the gap equation, and the minimum of the thermodynamic potential with 
respect to the quark condensates gives the "mass gap" equation, i.e. the equa- 
tion which controls the dependence of the dynamically generated quark masses. 

The energy eigenvalues cannot be computed analytically at all densities, 
except in two limiting situations. At low densities, where chiral symmetry is 
spontaneously broken, the gaps are zero. In this case, the quark dispersion 
relations are y^p^ + m*^ ± where m* = rrii — AGs {QiQi) are effective masses. 
The corresponding quark condensates are given by [95) 

<«">-ii>'*7?fc?^ (37) 

At high densities, the gaps are larger than the quark masses, hence the latter 
do not play significant roles and can be ignored. With progressively increasing 
density color-flavor locking becomes increasingly perfect, hence the name the 
CFL phase. In this phase, we can assume flavor symmetry and 8 of 9 quarks 
(3 colors times 3 flavors) have the dispersion relation \J {p ~ /i)^ + A^, while 
the remaining quark has the dispersion relation \J{j) — /-t)^ + 4A-^. In general, 
these properties are coupled so that both the masses and the gaps appear in the 
dispersion relations in a nontrivial fashion. 



Results of calculations based on the above model are displayed in Fig. 14 



using Eq. ( 35 1 . The dynamically generated quark masses are larger than the 
current quark masses at low density where chiral symmetry is spontaneously 
broken. In the CFL phase, the superconducting gaps form among all three 
combinations of unlike flavors, up-down, up-strange, and down-strange. In the 
2SC phase, the only pairing is between up and down quarks. This model exhibits 
a first-order phase transition between the gapped and ungapped phases, so 
the gaps do not continuously go to zero at low densities. The decrease in the 
gaps as a function of increasing density or large values of the quark chemical 
potential /j, is an artifact of the ultraviolet cutoff (a necessity imposed by the 
nonrenormalizable Lagrangian). HDET models show that the gaps increase 
with increasing [i. 



28 



400 



300 



T=0 



200 - 



100 



^ 
300 




400 500 
/X (MeV) 



100 



80 



60 



1 1 1 1 
- T=0 


1 1 1 1 1 
2SC 

K - 


/ 


7/ 




® / CFL 




us, 


ds / 




1 

1 

1 
1 



300 400 500 

At (MeV) 



Figure 14: Down and strange dynamical quark masses (left panel) and super- 
conducting gaps (right panel) as a function of density from [33]. The abscissa 
shows the quark chemical potential /x. The thin (bold) curves show values when 
color neutrality is not (is) enforced. Different results are obtained in the CFL 
and 2SC phases as shown. 



5.3 The Many Phases of Quark Matter 

There are many different possible pairing configurations in addition to the CFL 
and 2SC phases described above, including gapless phases [96], and color-spin 
locked pairing. Color superfluids also admit a new set of Goldstone bosons 
associated with flavor rotations of the pairing condensate which have a similar 
group structure to the pseudoscalar Goldstone bosons in QCD (tt, i^T, etc.). 
These bosons can condense [97] , forming a new phase of superconducting quark 
matter. The most common is the "CFL-K" phase which contains CFL quarks 
with a K meson condensate. All of these phases have their own associated 
exotic neutrino emissivities, including their own associated quark PBF neutrino 
cooling processes. A caricature phase diagram is shown in Fig. |15[ 

For densities near the deconfinement phase transition, the ground state of 
the quark superfluid may be similar to the Larkin-Ovchinnikov-Fulde-Ferrell [99| 
1100] (LOFF) pairing observed in condensed matter systems. LOFF pairing oc- 
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Figure 15: A schematic model QCD phase diagram, adapted from |98) . 



curs when two species participating in a pairing interaction have different Fermi 
momenta thus creating Cooper pairs with nonzero momentum. This pairing 
configuration breaks translational symmetry and encourages the formation a 
crystal. This creates a novel mechanism for pulsar glitches: the superfluid vor- 
tices pin to the crystalline part of the quark phase and during a glitch event the 
vortices move outwards by unpinning and repinning to the lattice |101j . This 
mechanism has not been either verified or ruled out by the data. 



5.4 Detecting Quark Matter 

The appearance of quark degrees of freedom often lowers the pressure at high 
densities, yielding neutron stars with smaller radii and smaller maximum masses 
compared to those in which quarks are absent. This is not guaranteed, however, 
as quark-quark interactions are sufficiently uncertain that quark matter can be 
nearly indistinguishable from matter consisting entirely of neutrons, protons, 
electrons and muons [102] . The presence (or absence) of quark matter will be 
difficult to determine from observations of neutron star structure alone. The 
neutron star mass-radius relation is dependent on the pressure of matter at a 
given energy density, but is insensitive to the nature of the particular degrees 
of freedom which provide that pressure. For this reason, it is natural to look to 
neutron star cooling to discern the composition of a neutron star's core. This 



topic will be considered in Sec. 7.7 
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6 Neutrino Emission Processes 



The thermal evolution ol neutron stars with ages ^ 10^ yrs is driven by neutrino 
emission. We wih here briefly describe the dominant processes; the interested 
reader can find a detailed description in 103 and an alternative point of view in 
|104j . Table [l] presents a short list of neutrino processes with estimates of their 
emissivities. Most noticeable is the clear distinction between processes involving 
5 degenerate fermions with a dependence, which are labeled as "slow" , and 
those with only 3 degenerate fermions with a dependence, which are several 
orders of magnitude more efficient and labeled as "fast" . The difference in the T 
dependence is important and is simply related to phase space arguments which 



are outlined in Sec. 6.5 The last subsection, § 6.6 describes the "PBF'"(Pair 
Breaking and Formation) process , which can provide a signal of the onset of 
a pairing phase transition that may be directly observable (and has likely been 
recently observed as described in Sec.|9]). 

6.1 The Direct Urea Process 

The simplest neutrino emitting processes are Eq. ^ and Eq. (|8| (see also Ta- 
ble [T]) , which collectively are generally referred to as the nucleon direct Urea 
("DU" or "DUrca") cycle. By the condition of /3-equilibrium, both reactions 
naturally satisfy energy conservation, but momentum conservation is more del- 
icate. Due to the high degree of degeneracy, all participating particles have 
momenta p{i) equal (within a small T <^ Tp correction) to their Fermi mo- 

1/3 

menta ppi'i)- As ppij-) oc and rip ^ <^ ri„ in neutron star matter, 

momentum conservation is not a priori guaranteed. In the absence of muons 
and hence with Up = rie, the "triangle rule" for momentum conservation requires 
that the proton concentration Xp > 1/9 ~ 11%, whereas at p ~ pnuc we have 
Xp ~ 4%. In the presence of muons, which appear just above Pnuc, the condition 
is stronger and one needs Xp larger than about 15% 105 . The proton fraction 
Xp grows with density (see Eq. (Ill), its growth being directly determined by 



the growth of the nuclear symmetry energy, so that the critical proton fraction 
for the DU process is likely reached at some supra- nuclear density |105j . For the 
APR EOS [16] that we will frequently use, the corresponding critical neutron 
star mass for the onset of the nucleon DU process is l.OTM©, but other EOSs 
can predict smaller critical densities and masses. 

It should also be noted that the direct Urea process, and for that matter, 
all the processes discussed in this section, can involve other leptons. Thus, for 
p ^ Pnuc, where /i^ appear, one also has 

n — >p + fj,^+D^, P + M^ — yn + i'f^. (38) 
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Table 1: A sample of neutrino emission processes. Tg 
of 10^ K and the i?'s are control factors to include 



is temperature T in units 
the suppressing effects of 



pairing (see See. 4.3) 
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6.2 The Modified Urea Process 



At densities below the threshold density for the nucleon DU process, where 
the DU process is forbidden at low temperatures, a variant of this process, 
the modified Urea ("MU" or "MUrca") process (see Table [l]) can operate, as 
advantage is taken of a neighboring nucleon in the medium |106j to conserve 
momentum. As it involves the participation of five degenerate particles, the 
MU process is much less efficient than the DU process. Unlike the nucleon DU 
process, which requires sufficient amount of protons, both branches of the MU 
process operate at any density when neutrons and protons are present. 

6.3 Bremsstrahlung 

Related to the MU processes is another class of processes, bremsstrahlung, made 
possible through neutral currents [f 07] . These differ from MU processes in that 
each reaction results in the production of a pair, and the pair can have any 
neutrino flavor. Bremsstrahlung reactions are less efficient, by about 2 orders of 
magnitude, than the MU processes, but may make important contributions in 
the case that the MU process is suppressed by pairing of neutrons or protons. 
Bremsstrahlung involving electron-ion scattering is also an important source of 
neutrino emission in neutron star crusts: 

e- + {A, Z) — >e- + {A,Z) + iy + iy, (39) 

where (A,Z) designates the participating ion. 

6.4 Exotic Matter: Hyperons, Deconfined Quarks, Meson 
Condensates 

In the presence of hyperons, DU processes which are obvious generalizations of 
the nucleon-only process, can also occur [59] and several are displayed in Table[l] 
When they appear, the A's initially have a density much smaller than that of 
the neutron and hence a smaller Fermi momentum. Consequently, momentum 
conservation in the A DU cycle is easily satisfied, requiring a A concentration 
XA ^ 3%. Notice that if the nucleon DU process is kinematically forbidden, the 
E~ DU process is also kinematically forbidden, whereas no-nucleon DU pro- 
cesses, of which one example in shown in Table [T] are possible. This particular 
no-nucleon DU process requires relatively low A and E~ threshold concentra- 
tions. Other examples involving S~, E'^, S"*", S~ and 'EP hyperons are given in 

m- 

In deconfined quark matter, DU processes involving all three flavors are 
possible, as indicated in Table [T] Rates for these processes have been calculated 
by Iwamoto |108| . 

Although not shown in Table [T] hyperons or quarks could also be involved in 
MU-like and bremsstrahlung processes (a quark MU process would involve an 
additional quark in the entrance and the exit channels, for example), but with 
greatly reduced rates compared to their corresponding DU processes and a 
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dependence. These processes are usually neglected since the DU processes are 
almost invariably allowed in the presence of hyperons or quarks. 

In the presence of a meson condensate, copious neutrino emission in the 
processes listed in Table [l] occurs [1091 I110| . As the meson condensate is a 
macroscopic object, there is no restriction arising from momentum conservation 
in these processes. 



6.5 Temperature Dependence of Neutrino Emission 

We turn now to briefly describe how the specific temperature dependence of 
the neutrino processes described above emerges. Consider first the simple case 
of the neutron /3-decay. The weak interaction is described by the Hamiltonian 
%! = {G F I \J2)B ^j,L^ , where Gp is Fermi's constant, and = '0e7'^(l — 75)V'i' 
and i?^ = V'p7Ai(C'yl ~ CAl^)i>n are the lepton and baryon weak currents, 
respectively. In the non relativistic approximation, one has B'^ = cos 6'j,^'],5'n 
and — ~ cos 6c Qa ^'p cr* '^n where 9c is the Cabibbo angle and gA the axial- 
vector coupling. Fermi's Golden rule gives us for the neutron decay rate 

- III - ^ l"'-!^ • 

i.e., a sum of (27r)^(5^(P/ — Pi) ■ \Mfi\'^ over the phase space of all final states / = 
{Pv,Pe,Pp)- The integration gives the well known result Wp = G|,cos^6'c(l + 
igjs^m^^c^wp / [211^), where wp ^ 1 takes into account small Coulomb correc- 
tions. This gives the neutron mean life, r„ ~ 15 minutes. Alternatively, a 
measurement of r„ determines Gp (modulo cos 0c and wp). 

The emissivity of the DU process (the Feynman diagram for this process 



is shown in Fig. 16 1 can be obtained by the same method, leading to 



(41) 

with an extra factor E^, for the neutrino energy and the phase space sum now 
includes the initial n. The fi terms, fi being the Fermi-Dirac distribution for 
particle i at temperature T, take into account: (1) the probability to have a n in 
the initial state, /„, and (2) the probabilities to have available states for the final 
e and p, denoted by (1 — /e) and (1 — /p), respectively. We do not introduce a 
Pauli blocking factor (1 — /jr) for the anti- neutrino as it is assumed to be able to 
freely leave the star (i.e., /jj = 0). When performing the phase space integrals, 
each degenerate fermion gives a factor T, as particles are restricted to be within 
a shell of thickness ksT of their respective Fermi surfaces. The anti-neutrino 
phase space gives a factor T^. The factors is ^ T and the delta function 
6'^{Pf — Pi) gives a factor from S{Ef — Ei). Altogether, we find that 

cx • T • T • T • ^ • (1)^ -T ^T^ , (42) 
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where the (1)^ factor eraphasizes that the squared matrix element \M 
T-independent. An exphcit expression for the neutrino emissivity for the DU 
process can be found in |105j . 




Direct Urea 



Modified Urea 



Bremsstralilung 



Figure 16: Feynman diagrams for the indicated neutrino emitting processes. 

Figure [16] shows a Feynman diagram for the MU process. There are two 
more such diagrams in which the weak interaction vertex is attached to one of 
the two incoming legs. In this case, the T-power counting gives 



(xT'^-T-T-T-T-T- 



1 

f 



{\f ■T = T^ 



(43) 



In this case, the involves two strong interaction vertices, connected by the 

wavy line in Fig. |16[ which is momentum independent and hence T-independent. 
The numerical difference ~ lO'^Tg between the MU and the DU rates, see 
Table [ij comes mostly from the extra phase space limitation cx from the 
two extra nucleons: as a dimensionless factor it is [UbT / EpY — 10~®T| with 
Ep ^ 100 MeV and fc^lO^ K ~ 0.1 MeV. Notice that in the MU case, the 
internal neutron is off-shell by an amount ~ /^e which does not introduce any 
extra T-dependence as we are working in the case Ep{e) ^ 100 MeV ^ T. 
Reference |106j contains the expression from which neutrino emissivity from the 
MU process can be calculated. 

Turning to the n—n bremsstrahlung process, one diagram is shown in Fig. |16| 
and there are three more diagrams with the weak interaction vertices attached 
to the other three external lines. The T-power counting now gives 



oc T'^ • T-^ • T • T • T • T ■ 



T 



(44) 



with two T^ factors for the neutrino pair. The factor (T^^)^ arises from the 
matrix element as the intermediate neutron is almost on-shell, with an energy 
deficit ~ T, and its propagator gives us a dependence for Mfi. A working 
expression for the bremsstrahlung process can be found in |106j . 
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6.6 The Cooper Pair Neutrino Process 

The formation of the fermonic pair condensate also triggers a new neutrino 
emission process [7TJ 11111 1112j which has been termed the "pair breaking and 
formation", or PBF, process |113j . Whenever any two fermions form a Cooper 
pair, the binding energy can be emitted a,s a, v — V pair. Under the right 
conditions, this PBF process can be the dominant cooHng agent in the evolution 
of a neutron star |114j . Such efficiency is due to the fact that the pairing 
phase transition is second order in nature. During the cooling of the star, the 
phase transition starts when the temperature T reaches when pairs begin to 
form, but thermal agitation will constantly induce the breaking of pairs with 
subsequent re-formation and possible neutrino pair emission. 



The emissivity of the PBF process (see the left panel of Fig. 17 for a Feynman 
diagram) can be written as 

^ //// 1:? (S^ ^ - ''■"'''•I' ■ 

(45) 

Under degenerate conditions, the expression above can be reduced to read as 

- ^SS^ (fcBT)^a,,i?, [A,(T)/T] 
= 3.51 X 10^1 ^ X rhf pFj Tl a/,, R, [A, (T)/T] (46) 

for a fermion / in a pairing state j = ^Sq or ^P2. The coefficients a/j depend 
on the type of fermion and on the vector and axial couplings Cy and Ca (see. 



e.g., [22]). The control functions Rj are plotted in the right panel of Fig. 17 
These functions encapsulate the effect that the PBF process turns on when T 
reaches Tc and practically turns off at T ^ 0.2 Tc when there is not enough 
thermal energy to break pairs. 

The PBF process has had an interesting history. It was first discovered by 
Flowers, Ruderman, and Sutherland [71j and, independently, by Voskresensky 
and Senatorov |lllj . It was, however, overlooked for 20 years, until implemented 
in a cooling calculation in |113j and its importance emphasized in '1141 . Then in 
2006, Leinson and Perez |115j showed that the previous computations of the PBF 
emissivity were incompatible with vector current (e.g. baryon number) conser- 
vation. Neutrino pair-production is mediated by the weak interaction, which can 
be decomposed in the traditional manner into vector and axial-vector parts. In 
pure neutron matter, the vector part of the PBF emissivity is suppressed be- 
cause of vector current conservation by a factor of order {pp/m*c)'^ — {vp/c)'^- 
This is equivalent to the simple observation that a one-component system of 
charges does not have a time- varying dipole moment |116j . The axial part of 
the PBF process is, however, unimpeded and dominates the emissivity. 
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Figure 17: Left panel: Feynman diagram for vV emission from the formation 
of a nn Cooper pair (pair breaking and formation, PBF, process). Right paneh 
control functions i?pBF for the PBF process. 



Temperature dependence of the PBF neutrino emissivity 



The temperature dependence of the PBF process (left panel of Fig. 17) can be 



ascertained from Eq. (45) according to the following T-power counting: 



^PBF . y3 . y . . _! . ^(A/T) • T = i?(A/T) , (47) 

where the two and the first T factors arise from the phase space integrations 
of the neutrino pair and the first participating nucleon, respectively. The fac- 
tor 1 results from the phase space integration of the second nucleon. As there 
are only two degenerate fermions in this process (in contrast to the Urea and 
bremsstrahlung processes that involve 3, 4, or 5 degenerate fermions), the mo- 
menta of the neutrino pair and the first nucleon are chosen the momentum of 
the second nucleon is fixed by the three-momentum conserving delta function. 
Thus, this second nucleon does not introduce any T dependence. The de- 
pendence arises from the energy conserving delta function. The last T factor is 
from the neutrino pair's energy, whereas the T and A dependence of the matrix 
element of the reaction are included in the function R{/S./T), which is shown in 
the right panel of Fig. [17] 

An alternative way of looking at the PBF process is simply as an interband 
transition of a nucleon |117j . Considering the particle spectrum in a paired state 
(the right panel of Fig. [s]), the lower branch (with e < ej? — A) corresponds to 
paired particles whereas the upper branch to excited ones, i.e., the "broken 
pair" leaves a hole in the lower branch. A transition of a particle from the 
upper branch to a hole in the lower branch corresponds to the formation of a 
Cooper pair. 
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Dominance of triplet-pairing 

In the non-relativistic limit for nucleons, the leading contribution from the axial- 
vector part is proportional to {vp/cf'- To this order, the control function R{z = 
A/T) receives a contribution from the axial- vector part which can be expressed 
as [m] 

2 r roc 6 J 

where x — vp{p — pp) /T and the quantity / = Ixx + lyy + Izz with 

= Y,{B\i>^^^^A){B\i>^au^A)* , (49) 
vv' 

where 5* is the second-quantized non-relativistic spinor wave function of the 
nucleons in superfluid matter (see |117) for its detailed structure in terms of 
the Bogoliubov transformation matrix elements Uariip) and Vo-^(p)), \A) is the 
initial state of the system and \B) its final state. The total spin states of the 
pair 77, rj' and a each take on values ±1 and Ui's are the Pauli spin matrices. 
The energy of a paring quasiparticle is given by i? = y^e^ + A^, where e = 

vf{p — Pf)- For singlet-state pairing, the momentum-dependent gap Ap is 
independent of p, so that the occupation probabilities Up and Vp associated 
with the matrix elements Ur,T)is>) and Vct^(p) depend only on p = |p|. In the 
case of singlet pairing, the Bogoliubov matrix elements satisfy the symmetry 
properties Vafi{—V) = ^^^(p) and Va^(p) = — Vsq(p), so that the diagonal 
elements of this 2x2 matrix are zero, and the non-diagonal elements are 
and —Vp^ respectively. These symmetry properties, together with the traceless 
property of the Pauli matrices, ensure that the quantity / = 0. Thus, to order 
(wf/c)^, the axial- vector part does not contribute in the spin-singlet channel 
rendering the triplet pairing channel, which does not vanish, to be the sole 
contribution to the PEP process. 

Time history of the PBF process 

The Cooper pair neutrino process operates at different times in a neutron star's 
cooling history according to the time during which the local temperature is 
nearly equal to the critical temperature of any superfluid gap. In neutron stars 
consisting of neutrons, protons, and electrons, there are three relevant super- 
fluid gaps: singlet neutron superfluidity, singlet proton superfluidity, and triplet 
neutron superfluidity at high densities when the singlet channel of the neutron- 
neutron interaction becomes repulsive above the saturation density. In neutron 
stars which contain exotic matter in their interiors, each additional superfluid 
fermion potentially opens up new Cooper pair cooling processes. If neutron stars 
contain deconfined quark matter in their cores, then pairing between quark fla- 
vors creates new Cooper pair neutrino processes which involve pairing between 
unlike fermions |118) . 
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7 Cooling of Neutron Stars 



The study of neutron star cooling is a Sherlock Holmes investigation, following 
the tracks of energy. At its birth, some 300 B (1 Bethe = 10^^ ergs) of grav- 
itational energy are converted largely into thermal energy. About 98% of it is 
emitted in neutrinos during the first minute, the proto-neutron star phase, 1% 
is transferred to the supernova ejecta (with 1% of this 1% powering the light 
show), and the remainder is left in thermal energy of the new-born neutron star, 
i.e., the star produced during the proto-neutron star phase. Following the tracks 
of energy, the subsequent evolution of the neutron star can, in a simplified way, 
be described by an energy balance equation 

5!§i=ft.f ^-L„-L, + i,, (50) 

where Eth is the star's total thermal energy, Cy its specific heat, and Lj and 
Li, its photon and neutrino luminosities, respectively. The term H, for "heat- 
ing" , represents possible dissipative processes, such as friction from differential 
rotation or magnetic field decay. In this simplified equation it is assumed that 
the star's interior is isothermal with temperature T, a state reached within a 



few decades after birth in the core-collapse supernova (see Sec. 7.4). A more de 



tailed study would include general relativistic effects and consider a local energy 



balance equation for each layer in the star, instead of the global one of Eq. ( 50 ), 
complemented by a heat transport equation, in order to follow the evolution 
of the temperature profile in the stellar interior (see, e.g., [HT] and references 
therein) . 

After the proto-neutron star phase, matter is highly degenerate within most 
of the star, except the outermost, lowest density, layers. As a consequence, 
the gross structure of the star does not evolve with time and is determined, 
once and for all, by solving the Tolman-Oppenheimer-Volkoff equations [13] 
of hydrostatic equilibrium. An equation of state is required to not only solve 
these equations, which determine the mass and radius of the star, but also 
to evaluate the internal chemical composition of each species of nucleus and 
particle as well as their effective masses, chemical potentials, specific heats, etc. 
A complete cooling model requires, moreover, inclusion of neutrino and surface 
photon emissions as well as a description of the pairing properties of matter, 
i.e., the pairing gaps for each fermonic species, together with their respective 
density dependences. 



Within the isothermal approximation of Eq. (50), the three major ingre- 
dients needed for the study are Cy, L-y, and L^. Neutrino emission processes 
were described in Sec.|6]and the specific heat and photon emission are briefly de- 
scribed below. We continue this section by describing simple analytical solutions 



of Eq. ( 50 ) and displaying the results of representative numerical simulations of 



the complete set of general relativistic evolutionary equations. 
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Figure 18: Cumulative specific heats of e, /j,, p, and n as a function of stellar 
volume within the core of a 1.4 Mq star built using the APR EOS at T — 10^ 
K. Nucleons are assumed to be unpaired which implies Cv oc T. No hyperons or 
quarks are permitted by the EOS. This figure is adapted from 



7.1 Specific Heat 

The dominant contributions to the specific heat Cv come from the core, which 
makes up more than 90% of the total volume and 98% of the mass. Its con- 
stituents are quantum liquids of leptons, baryons, mesons, and, possibly, decon- 
fined quarks at the highest densities. Hence, one has 

Cv = ^Cv,i with Cv,i = Cyidv, (51) 



where Cv,i is the specific heat per unit volume of constituent i (i = e, fj,, n, p, 
hyperons, quarks), but those of meson condensates is usually neglected.. For 
normal (i.e., unpaired) degenerate fermions, the standard Fermi liquid result 

mi 

Cv.-iV(0)^4r with iV(0) = "^ (52) 



can be used, where m* is the fermion's effective mass. In Fig. 18 the various 
contributions to Cy are illustrated. 

When baryons, and quarks, become paired, as briefiy described in Sec. |4.3[ 
their contribution to Cy is strongly suppressed at temperatures T <^ {T^ 
being the corresponding critical temperature). Extensive baryon, and quark, 
pairing can thus significantly reduce Cy, but by no more than a factor of or- 
der ten because leptons do not pair. The crustal specific heat is, in principle, 
dominated by neutrons in the inner crust but, as these are certainly extensively 
paired, only the nuclear lattice and electrons contribute in practice. 
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7.2 Photon Thermal Luminosity and the Envelope 

The photon thermal luminosity is commonly expressed through the effective 
temperature defined by 

= AirR^asBTt , (53) 

where (tsb is the Stefan-Boltzmann constant. Thermal photons from the neu- 
tron star surface are effectively emitted at the photosphere, which is usually in 
an atmosphere, but could be located on the solid surface if a very strong mag- 
netic field exists [T7]. The atmosphere, which is only a few centimeters thick, 
contains a temperature gradient; gives an estimate of its average tempera- 
ture. The opacity in the atmosphere receives a strong contribution from free- free 
scattering that has a strong E^^) energy dependence. As a result photons of 
increasing energy escape from deeper and hotter layers and the emitted thermal 
spectrum shows an excess of emission in its Wien's tail compared to a black- 
body of the same temperature T = Tg |120j . In the presence of heavy elements, 
"metals" in astronomical parlance, which may be not fully ionized, absorption 
lines increase the bound-free opacity contributions, and push the Wien's tail of 
the observable spectrum closer to the blackbody one. The presence of a strong 
magnetic field also alters the opacity in such way as to mimic a blackbody with 
the same chemical composition and Tg. There were great expectations that, 
with the improved spectral capabilities of Chandra and Newton observatories, 
many absorption lines would be observed and allow the determination of the 
gravitational redshifts and chemical composition of isolated neutron star sur- 
faces. This expectation has, unfortunately, not been fulfilled; only in a very few 
cases have lines been detected, and their interpretation is controversial. 



Observationally, and Te Sltq red-shifted and Eq. (531 is rewritten as 

L^=47ri?^asB(rr)\ (54) 

where = e^'^L^, = e'^T^, and = e-'t'R. Here e'^ = 1 + z, with 
z being the redshift, and e^'^ is the goo coefficient of the Schwarzschild metric. 



Notice that i?oo bas the physical interpretation of being the star's radius corre- 
sponding to its circumference divided by 2tt, and would be the radius one would 
measure trigonometrically, if that were possible [121]. 

and are the observational quantitities that are compared with the- 
oretical cooling models. In principle, both are independently observable: 
is deduced from a fit of the observed spectrum while L'^ is deduced from the 
observed total fluxj^ and knowledge of the distance D, via = AirD^F. The 
star's distance can be deduced either from the radio signal dispersion measure, 
if it is a radio pulsar, or from the distance of an associated supernova remnant. 



°Tho flux must, however, be corrected for intersteUar absorption. 
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if any. Then Eq. ( 53 1 or ( 54 1 provides a consistency check: the inferred radius 
R should be of the order of 10 - 15 km. Given the lack of determination of the 
atmospheric composition from spectral lines, this consistency check is generally 
the only criterion to decide on the reliability of a measurement from an at- 
mosphere model spectral fit (besides the obvious requirement that the model 
must give a good fit to the data, i.e. a — !)• 

In a detailed cooling calculation, the time evolution of the temperature pro- 
file in the star is followed. However, the uppermost layers have a thermal time- 
scale much shorter than the interior of the star and are practically always in a 
steady state. It is, hence, common to treat these layers separately as an enve- 
lope. Encompassing a density range from pb at its bottom (typically pf, — 10^° 
g cm"-^) up to pe at the photosphere (pe ^ 1 g cm""^), and a temperature range 
from Th to Tg, the envelope is about one hundred meters thick. Due to the 
high thermal conductivity of degenerate matter, stars older than a few decades 
have an almost uniform internal temperature, except within the envelope which 
acts as a thermal blanket insulating the hot interior from the colder surface. A 
simple relationship between and Tg can be formulated |122j : 



^6 f 



0.5+a 



with a <C 1. The precise Te — Tj, relationship depends on the chemical com- 
position of the envelope. The presence of light elements like H, He, C, or O, 
which have large thermal conductivities, leads to a larger for the same Tf, 
relative to the case of a heavy element, such as iron, envelope. Light elements 
are not expected to survive densities larger than ~ 10^° g cm~^ due to pycnonu- 
clear reactions. Thus, the maximum possible mass in light elements amounts to 
AMLight — 1O~^M0, which is enough to raise Tg by a factor of two. Magnetic 
fields also alter the Te — relationship, but to a lesser extent (see, e.g., |123) 
for more details) unless they are super-strong as in the case of magnetars, i.e. 
B, ~ 10^^ G. 



7.3 Analytical Solutions 



As the essential ingredients entering Eq. ( 50 ) can all be approximated by power- 
law functions, one can obtain simple and illustrative analytical solutions (see 
also |124p . We adopt the notation 

Cy = CgTg , = TVgTl , and = SgT^+^'' , (57) 

where Tg = T/(10^ K) refers to the isothermal temperature Ti, in the star's 
interior. As written, L^, considers slow neutrino emission involving five degener- 
ate fermions from the modified Urea and the similar bremsstrahlung processes, 
summarized in Table [I] The photon luminosity is obtained from Eq. ( 53 1 



using the simple expression in Eq. (56 1. We will ignore redshift. Typical pa- 
rameter values are Cg ~ 10^^ erg K~, iVg ~ 10^° erg s~^ and Sg ~ lO'^'^ erg 
s~^ (see Table 3 in [123j for more details). In young stars, neutrinos dominate 
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the energy losses (in the so-called neutrino cooling era), and photons take over 
after about 10^ years (the photon cooling era). 

Neutrino cooling era: In this case can be neglected in Eq. (50), so that 



rp6 Ti6 

^9 ^0,9 



(^)''' (When r« To) 



(58) 



with a MU cooling timescale tmu = IQ^Cg/GN^ ~ 1 year when the star reaches 
the asymptotic stage (T <C Tq, Tq being the initial temperature at time t = 0). 
The observed slope of the cooling track during the asymptotic stage is 



dlnTp 



dint 



1 a 
12 6 



(59) 



using the core-envelope relation Eq. ( 56 1 . This slope is not sensitive to the 
core-envelope relation. 

Photon cooling era: In this era, is negligible compared to L-y. Since 
|a| <C 1, one finds 



t^ti+T.^ [In Tq - In Ti^g] Tg = Ti, 



{ti-t)/r-, 



(60) 



where Ti is T at time t = ti and Tj = lO^Cg/Sg ~ 3 x 10^ years, the photon 
cooling timescale. The observed slope of the cooling track is, when T <C Ti, 



dlnTp 



d\nt 



2t^ 



(61) 



This slope becomes steeper with time. 



7.4 Numerical Simulations and the Effects of Pairing 

Numerical simulations of a cooling neutron star use an evolution code in which 
the energy balance and energy transport equations, in their fully general rel- 
ativistic forms, are solved, usually assuming spherical symmetry and with a 
numerical radial grid of several hundred zones |^ Compared to the analytic so- 



lutions discussed in Sec. 7.3 detailed conductivities, neutrino cmissivities, and 
possible pairing are taken into account. 

A set of cooling curves that illustrate the difference between cooling driven 
by the modified Urea and the direct Urea processes is presented in Fig. [T9j 
Cooling curves of eight different stars of increasing mass are shown, using an 
equation of state model from [125j . which allows the DU process at densities 
above 1.35 x 10^^ g cm~'^, i.e., above a critical neutron star mass of 1.35 Mq. 
Notice that the equation of state used has parameters specifically adjusted to 
obtain a critical mass of 1.35 Mq which falls within the expected range of 
isolated neutron star masses; other equations of state can result in very different 



^Such a code, NSCool, is available at: http://www.astroscu.unain.mx/neutrones/NSCool/ 
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Log Age (yrs) 



Figure 19: Cooling curves illustrating the difFcrcnee between slow cooling driven 
by the modified Urea process, for masses below 1.35Mq, and fast cooling from 
the direct Urea process in more massive stars [57] . 



critical masses. The difference arising from slow and fast neutrino processes is 
clear. The various cooling phases, A to D are discussed below. 

To further illustrate some of the possible cooling behaviors of neutron stars, 
and the effects of pairing, we show in Fig. 20 simulations based on a IAMq 
star built with the APR EOS [TB] and a heavy-element envelope. The "slow 
cooling" models include, in the core, the slow neutrino processes of Table [T] 
and the PBF process only. For the "fast cooling" models, a fast process with 
emissivity e^']^j = lO' • erg cm s , with q — 25, 26, and 27, was added at 
p > 3pnuc- These q values simulate neutrino emission from a kaon condensate, 
a pion condensate, or a direct Urea, respectively. These models, being based 
on the same EOS, are not self-consistent, but they have the advantage that 
the only differences among them is the presence or absence of the fast cooling 
process with e^^^^ and the presence or absence of nucleon pairing. The models 
with pairing include the neutron ^Sq gap "SFB" of Fig. 10 in the inner crust, 
^Sq proton gap model 



the 
3p-F 



T" of Fig. 

2 gap "b" of the right panel of Fig 



11 and the phenomenological neutron 



12 



The distinctive phases of evolution are labelled "A" , "B" , "C" , and "D" on 
the cooling curves in Fig. [19] and above the cooling curves in the left panel of 
Fig. [20] Phases A and B are determined by the evolution of the crust while C 
and D reflect the evolution of the core. We describe these four phases in more 
detail according to Fig. 20 
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Phase A: The effective surface temperature Tg here is determined by the evo- 
lution of the outer crust only. At such early stages, the temperature profile in 
the outer crust is independent of the rest of the star and, as a result, all models 
have the same Tg. 



Phase B: The age of the star during this phase is similar to the thermal relax- 
ation timescale of the crust, the heat flow in which controls the evolution of Tg. 
The evolution of the temperature profile for the fast cooling model with q = 26 
in the absence of pairing (marked as "Normal" in the left panel of Fig. 20 1 is 
depicted in the right panel of this figure |20j Very early in the evolution, a cold 
"pit" develops in the core where fast neutrino emission is occurring. During the 
first 30 years, heat flows from the outer core into this pit until the core becomes 
isothermal. Afterwards, heat from the crust rapidly flows into the cold core and 
the surface temperature T^, drops rapidly. Well before 300 years, during phase 
C, the stellar interior becomes isothermal and it is only within the shallow en- 
velope, not shown in this figure to preserve clarity, that a temperature gradient 
is still present. 

Notice that models of with pairing have shorter crust relaxation times due 
to the strong reduction of the neutron specific heat in the inner crust by "'^Sq 
neutron superfluidity there. 



Phase C: This is the "neutrino cooling phase" in which the star's evolution 
is driven by neutrino emission from the core: ^ L^. The difference be- 
tween "slow" and "fast" neutrino emission, with or without core superfluidity, 
is clearly seen. 

Quite noticeable is the effect of pairing-induced suppression of the neutrino 
emissivity in the fast cooling models. Once T drops below Tc, which happens 
only a few seconds or minutes after the star's birth, neutrino cooling is quenched. 
Which fast cooling process occurs is much less important in the presence of 
pairing than in its absence. It takes half a minute, ii q = 27, or half an hour, if 
g = 25, for the "pit's" temperature to fall below Tc, which does not matter when 
looking at the star thousands of years later. The evolution is more dependent 
on Tc than on q. 

In the case of the slow cooling models, the effects of pairing are more sub- 
tle than that in the fast cooling models, if one ignores the artificial case with 
the PBF processed turned off. The burst of neutrino emission occurring when 
T ~ Tci from the PBF process, induces an additional, if temporary, rapid cool- 
ing episode. The impact of the PBF process, however, depends on the size of the 
neutron ■^P-F2 gap. If the gap is large enough, the PBF cooling occurs during 
stage B and is hidden; if the gap is small enough, the PBF cooling occurs during 
stage D and is again hidden. Only intermediate size gaps reveal the presence of 
the PBF process. The effect of this gap on the evolution is considered in more 
detail in Sec. [s] The effect of the proton ^So gap is more subtle still, and its 
effects are considered in detail in Sec. M 
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A B C D 




Figure 20: Left panel: cooling curves for various illustrative cooling scenarios. 
Right panel: temperature profile evolution for the fast cooling model with q = 
26. The numbers on the curves give the age of the star, in years. See text for 
description. Figure taken from |126) . 



Phase D: At late times, has decreased significantly due to its strong T 
dependence and photon emission, L^, which is less T-dependent, now drives the 
evolution. This is reflected by the larger slopes of the cooling curves. During 
this "photon cooling era" , models with pairing cool faster due to the reduction 
of the specific heat from superfluidity. 



The effect of a very strong magnetic field 

In this contribution, we have neglected effects of the presence of a magnetic 
field. A strong magnetic field can alter the cooling of a neutron star in two 
ways. In the envelope and the crust, heat is transported by electrons. A surface 
magnetic field of strength ^ 10^^ G, a typical value for the majority of pulsars, 
is sufficient to induce anisotropy in the thermal conductivity of the envelope 
resulting in a non-uniform surface temperature [I21J that manifests itself as a 
modulation of the thermal flux with the pulsar's rotational period. As mentioned 



in Sec. 7.2 this effect only slightly alters the — Tg relationship of Eq. (561 
|127j . However, a strong magnetic field deep in the crust |128 will have larger 
effects. For example, a ^ lO^^'G toroidal field within the crust can act as an 
efficient insulator, rendering most of the star's surface very cold, but having two 
hot spots on the symmetry axis of the torus [1291 1130j . This results in peculiar 
cooling trajectories jl31j . The second effect of a magnetic field is that a slowly 



decaying field can act as a source of energy (i.e., the "H" term in Eq. (50)) 
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which can keep old neutron stars warm. If the magnetic energy reservoir is large 
enough, and small scale magnetic structures exist that can decay rapidljF^l 




aided by the Hall drift, the thermal evolution is significantly altered |132) . 

7.5 Observations of Cooling Neutron Stars 

It is useful at this juncture to compare these general behaviors with observations. 
We summarize the observational information relevant to neutron star cooling in 
Fig. 21 The data is separated into three subsets of stars. In the upper panel of 
Fig. 21 are presented 13 stars for which a thermal spectrum, in the soft X-ray 
band (0.1 - 3.0 keV) is clearly detected. In the lower panel, we show data from 4 
pulsars, labeled (A) through (D), that arc seen in the X-ray band with a power- 
law spectrum, but whose detected emission is of magnetospheric origin. Since 
the surface thermal emission from these 4 stars is undetected, being covered by 
the magnetospheric emission, only upper limits on their effective temperatures 
could be inferred. Finally, the lower panel of this figure contains 6 upper limits 
resulting from the absence of detection of any emission from compact objects 
in 6 gravitational collapse supernova remnants, labeled (a) through (f), from a 
search by Kaplan and collaborators [133] . Since no compact object has been 
detected in these supernova remnants, some of them may contain an isolated 
black hole, but that is unlikely to be the case for all 6 of them. The observations 
and estimates of ages and temperatures, or upper limits, are detailed in [5P and 
updated in [52]. 

Estimates of surface temperatures (redshifted to infinity) require atmospheric 
modeling in which three factors are involved: (1) the composition of the atmo- 
sphere (i.e., H, He or heavy element); (2) the column density of X-ray absorbing 
material (mostly H) between the star and the Earth; and (3) the surface grav- 
itational redshift. The column density is important because the bulk of the 
emitted flux is absorbed before reaching the Earth. The redshift does not af- 
fect blackbody models, but can influence heavy-element atmosphere models. In 
most cases, this quantity is not optimized, but set to the canonical value for 
M = IAMq and i? = 10 km stars. 

Narrow spectral lines are not observed in any of these sources, so their atmo- 
spheric compositions are unknown. However, some information can be deduced 
from the shape of the spectral distribution as heavy-element atmospheres closely 
resemble blackbodies. Fitting the flux and temperature of a source to a model 
yields the neutron star radius, if the distance is known. In some cases, clothing 
a star with a light-element atmosphere results in a predicted radius much larger 
than the canonical value, and one can infer the presence of a heavy-element 
atmosphere. Chang & Bildsten |134j have noted from such radius arguments 
that there exists a trend for stars younger than 10^ years to be better fit by 
light-element atmospheres and stars older than 10^ years to be better fit by 
heavy-element atmospheres. The possible evolution of stars leading to this trend 
is discussed further in (5lj . 

^''The Ohmic decay time, e.g., is oc , I being the length-scale of the field structure. 
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Figure 21: The present data set of cooling neutron stars. See text for presenta- 
tion. 
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Before embarking on the presentation of detailled cooling scenarios and com- 
paring them with data, it is worthwhile to discuss one general feature here. The 
four oldest stars, numbered 10 to 13 in Fig. [2l] will apprear to be much warmer 
than most theoretical predictions. These stars are prime candidates for the 



occurence of some late "heating", i.e. the "i?" term in Eq. (50). Models can 
be easily adapted to incorporate such heating (which has nothing to do with 
neutrino emission or pairing) that may become important at these times 



We will therefore concentrate on younger objects in the comparison with data. 

7.6 An Example: Fast Cooling with Pairing 

The thermal evolution of a neutron star undergoing fast neutrino emission was 
illustrated in Fig. [20] which showed that such a star, after the initial crust 
relaxation phase, has a very low Tg unless the neutrino emission is suppressed 
by pairing. In the presence of a gap that quenches fast neutrino emission, the 
resulting Tg is more sensitive to the size of the quenching gap than to the specific 
fast neutrino process in action. 

In Fig. [22j we compare one specific scenario with the observational data 
described in the previous Sec.|7.5| This scenario employs the same EOS as in the 



models of Fig. 19 and is more realistic than what was described in Sec. 7.4 insofar 
as it takes into account uncertainties concerning the chemical composition of 
the envelope. If the neutron star has an envelope containing the maximum 
possible amount of light elements, (~ 10~^ Mq), its is raised by a factor 
of two compared to the case when its envelope contains only heavy elements. 
Results for these two extreme scenarios of envelope composition are separately 



shown in the two panels of Fig. 20 In each panel, the first noticeable feature is 
the mass dependence of results when the mass exceeds 1.35M0 and the direct 
Urea process is activated. This is partially due to the increasing size of the fast 
neutrino emitting "pit" as the star's mass increases, but more dominantly due 
to the decrease in the neutron ■^P-F2 gap with increasing density. More massive 
neutron stars, at least in this model, will have the direct Urea process acting in 
their inner cores for a longer time until their central temperatures drop below 
the corresponding central values of Tg. Notice that in a model where the '^P-F2 
gap keeps growing with density as, in the "3bf ' model of Fig. [12] this mass 
dependence essentially disappears. 



When comparing the cooling trajectories of Fig. 22 with the displayed data 
we find that, with the exception of some of the oldest stars, all observed tem- 
peratures can be reasonably well fit. This result implies that different neutron 
stars have different masses, and some stars require an envelope containing light 
elements. This conclusion is common to all fast cooling scenarios that involve 
either the simple direct Urea process with nucleons or any direct Urea process 
with hyperons or quarks, as well as scenarios with charged mesons condensates. 
The successful recipe consists of having some fast neutrino process allowed in 
massive enough neutron stars. Then, in order to prevent intermediate mass 
stars from cooling too much, it is necessary to have a gap, or several gaps, large 
enough so that all fast neutrino processes are quenched at some point. 
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Figure 22: Comparision of a cooling scenario with data. The dense matter EOS 
used to build the stars is a PAL [125j version [57] that has only nucleons and 
leptons, and allows for the direct Urea process to occur in stars of mass larger 
than = 1.35M0. Cooling trajectories are shown for seven different masses 
and two different envelope chemical compositions. The neutron '^P-F2 gap from 
is implemented while the neutron and proton ^Sq gaps are 

respectively. 



model "b" of Fig. 



12 



from models "SFB" of Fig. 10 and "T" from Fig. 11 
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7.7 Can Quark Matter Be Detected from Cooling Obser- 
vations? 



The previous subsection presented one specific fast cooling scenario, wliile Sec. 
|7.4| only described very general trends of fast cooling scenario with pairing. It 
is interesting, however, to briefly consider if observations could determine not 
only whether or not fast cooling has occurred but also which fast cooling mech- 
anism operates. The first study of cooling of hybrid stars with superconducting 
quarks |135j showed, in fact, that it will be difficult to establish that neutron 
stars contain deconfined quark matter based solely on the cooling data of iso- 
lated neutron stars. When quark superconductivity is present, for example in 
the 2SC phase, neutrino emission is suppressed by a factor of e^^^^ for each 
flavor and color quark which is paired. Cooling in the CFL phase is special be- 
cause of the Goldstone bosons, which provide new cooling processes, and other 
pairing configurations also have their own cooling behaviors. Quark gaps that 
are typically larger than baryonic gaps therefore imply that hybrid star cooling 
is driven by other degrees of freedom present at lower densities. On the other 
hand, small quark gaps, or the absence of quark matter entirely, imply similar 
cooling behaviors that are easily reproduced by models with varying masses. 



The above conclusions are demonstrated in Fig. 23 in different ways in 
the two panels, using two pairs of EOSs. In the left panel, cooling curves 
with a fixed neutron star mass and employing two EOSs, one containing only 
nucleons (np) and the other with nucleons with quark matter (npQ), are shown. 
These simulations utilize a range of nucleon and quark pairing gaps. In the 
right panel, cooling curves of neutron stars of various masses and employing 
two EOSs, one with nucleons and hyperons (npH) and the other with quark 
matter added (npHQ), are compared. This second set of simulations utilizes 
fixed nucleon/hyperon/quark pairing gaps. 

We refer the reader to [135] for details, but the figure illustrates that, given 
our poor knowledge of the size and density extent of the various pairing gaps, 
nearly indistinguishable cooling curves can be obtained with very different dense 
matter models, i.e with or without hyperons and/or with or without quark 
matter. 



7.8 Can Superfluid Gaps Be Measured, or Predicted? 

We are led, by the discussions in Sec. |7.6| and Sec. |7.7[ to an intriguing result: 
Any equation of state that permits fast neutrino emission is compatible with the 
presently available cooling data on isolated neutron stars IF pairing is present 
with gap(s) of the appropriate size(s)^^ This result, unfortunately, is ambigu- 
ous: it means that theoretical uncertainties in the description of dense matter 
make it difficult to determine the composition of neutron star cores from isolated 

^^The ultimate theory of neutron star matter will yield the equation of state and the super- 
fluid gaps. However, at present, this is far from being the case since the gaps, as described in 
Sec. [4] are very sensitive to many intricate Fermi surface processes to which the equation of 
state, which is a bulk property, is quite insensitive. 
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Log Age [yrs] Log Age [yrs] 

Figure 23: Left panel: cooling of a IAMq neutron star composed only of 
nucleons (continuous lines) or nucleons and quark matter (dashed and dotted 
curves), with various assumptions on the pairing gaps. Right panel: mass depen- 
dent cooling histories for two dense matter scenarios: "npH" matter contains 
nucleons and hyperons; "npHQ" matter also includes deconfined quarks. In 
these two scenarios, pairing gaps are assumed fixed. Figures from [135) . 

neutron star cooling observations alone. 

One could, however, take an optimistic approach and reverse the argument. 
As the cooling of neutron stars is very sensitive to the presence, and the sizes, of 
pairing gaps, we can apply this fact to attempt to measure the gap(s) [ST]- This 
is in line with the idea of using neutron stars as extra-terrestrial laboratories to 
study dense matter. However, as massive neutron stars may contain hyperons, 
meson condensates or deconfined quarks, the question arises: which gap would 
we be measuring? |135| . This situation, as demonstrated in the maximally 
complicated models of Fig. |23[ may appear to render further studies of neutron 
star cooling a waste of time. 

However, our discussion so far has assumed the presence of "fast" cooling. 
Is fast cooling, in fact, actually necessary? This question motivates the minimal 
cooling scenario, to which we turn in Sec. |8j which considers only models in 
which enhanced cooling processes do not occur. 

8 Minimal Cooling and Superfluidity 

The comparison in the previous section of theoretical cooling models with avail- 
able data suggests the presence of extensive pairing in the cores of neutron stars. 
However, the plethora of possible scenarios makes it difficult to go beyond this 
generic conclusion. And, in spite of the fact that many theoretical models of 
dense matter predict the presence of some form of "exotica", one should still 
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ask the question "do we need them?" and, if yes, how strong is the evidence for 
them? To address these questions, the minimal cooling paradigm was developed 
in [5T1[52]- In this scenario, all possible fast neutrino emission processes (from 
direct Urea processes involving baryons or quarks, and from any form of exot- 
ica) are excluded a priori. Superfluid effects along with the PBF process are, 
however, included in the slow neutrino emission processes, such as the modified 
Urea processes involving neutrons and protons. This is a very restrictive sce- 
nario that minimizes the number of degrees of freedom, but fully incorporates 
uncertainties associated with the equation of state, envelope composition and its 
mass, etc. A detailed presentation of this scenario can be found in [511 and a 
variant of it was developed by the loffe group in [1361 1137] . Extensive studies of 
these two groups have pinned down two major sources of uncertainties: (i) our 
present lack of knowledge on the chemical composition of the envelope, and (ii) 
the size and extent of the neutron •^P-F2 gap. It turns out that, in the absence of 
any "fast process" , neutrino emissivity resulting from the PBF process involving 
^P-F2 Cooper pairs is a major factor in the thermal evolution of a neutron star 
and its effect is strongly dependent on the size of the neutron '^P-F2 gap. We 
describe in some detail effects of the PBF process in the next subsection before 
turning to a comparison of the predictions of the minimal cooling paradigm with 
available data. 



8.1 Effects of the PBF Induced Neutrino Process 

The PBF process is distinctive in the sense that, in any layer in the star's 
interior, it turns on when the temperature T reaches the corresponding Tc of 
the ambient density. Then the pairing phase transition is triggered and while 
T is not too much lower than Tc, there is constant formation and breaking 
of Cooper pairs induced by thermal excitation. However, when T falls below 
^ 0.2 Tc, there is not enough thermal energy to break pairs and the process 



shuts off (sec its control function displayed in the right panel of Fig. 17). As 



described in Sec. 6.6 the neutrino emissivity of the PBF process is significant 
only in the case of anisotropic pairing, as in the case of a '^P-F2 gap. 

The schematic Fig. [24] illustrates the effect of the pairing phase transition 
on the total neutrino luminosity of the star, L,^, depicted in the upper part of 
each panel as a function of the core temperature T. The long-dashed curves 
labelled "MUrca" show the modified Urea luminosity, cx in the absence 
of neutron ^P-F2 pairing while the short-dashed lines labelled "PBF" show the 



maximum possible PBF luminosity, L^^^ oc T^, i.e. from Eq. (46) with the 
control function Rj set to 1. The lower part of each panel shows a neutron 
^P-F2 Tc versus radius curve (displayed with its T axis horizontal to coincide 
with the upper part's T axis). Assuming the temperature T to be nearly uniform 
in space in the entire core, the core T profile is just a vertical straight line that 
moves from the right to the left as the star cools down. Each panel considers a 
different type of Tc curve and in its upper part, the thick continuous line shows 
the actual evolution of the total in the presence of pairing. As the star cools, 
and while T > T^^, the total Li, is dominated by the MUrca process with small 



53 



max 

cn ' 



contributions from other much less efficient processes. When T reaches 
there is a thick layer in which the phase transition starts and where neutrino 
emission from the PBF process is triggered which increases L^. As T decreases, 
but is still larger then T™'', there will be one or two layers in the star where 
T is only slightly lower than the corresponding value of in that layer(s) and 
where the PBF process is efficiently acting. (At these times, neutrino emission 
is strongly suppressed in the layers where T is much smaller than Tc.) Finally, 
when T drops well below T)!^™, neutrino emission by the dominant processes is 
suppressed everywhere in the core and drops rapidly to reflect the much less 
efficient neutrino emissivity from the crust. The higher the value of the 
earlier in the star's history will the enhanced cooling from the PBF process be 
triggered, while the lower the value of T™™ the longer this enhanced cooling will 
last, as illustrated by the difference between the left and right panels of Fig. 24 
This description will be pursued with a simple analytical model in Sec. |9.1|and 





Figure 24: Schematic illustration of the PBF process induced neutrino luminos- 
ity as controlled by the shape of the Tc curve. See text for a detailed description. 



From the above discussion, we learn that the effect of the PBF process 
depends strongly on the size and density dependence of the neutron '^P-F2 gap. 
To explore the possible impact of these two physical ingredients, we show in 



Fig. 25 two sets of cooling trajectories resulting from different assumptions. 
For contrast, both panels show a curve with no pairing, i.e., the same as the 
slow cooling model marked "Normal" in the left panel of Fig. [20] The various 
cooling curves are labelled by T™^'^, the maximum value of reached in the 
stellar core. In the left panel, large values of corresponding to the three 

phenomenological gaps "a", "b", and "c" of the right panel of Fig. [121 are 
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employed. All three gaps lead to the same result during the photon cooling 
era {t > 10^ yrs), but case "a" strongly differs from the other two cases during 
the neutrino cooling era. All three models have large enough gaps and T^'^^ 
that the enhanced cooling from the PBF process started well before the end of 
the crust relaxation phase. The gaps and T™" of models "b" and "c" are so 
large that they induce a very strong suppression of at early times and hence 
high values of Tg during the neutrino cooling phase. In contrast, the gap "a" 
has a moderate T™'^'^ and also a very small T™" which guarantees that at ages 
10^ — 10^ yrs there is alway a significant layer going through the pairing phase 
transition in which L^^^ is large which results in a colder star. 

As small gaps have the strongest effect in inducing some enhanced cooling of 



the neutron star, we explore in the right panel of Fig. 25 the effect of increasingly 
small neutron '^P-F2 gaps by scaling down the model "a" gap by a factor 0.6, 
0.4, and 0.2, respectively. A new feature now emerges: with decreasing values 
of T™'^^, the pairing phase transition is initiated at progressively later stages. 
This feature is signaled by cooling trajectories (dotted curves) departing from 
the T^'^^ = trajectory. This late onset of pairing then manifests itself as a 
sudden rapid cooling of the star (due to the sudden increase of L^, from 
to L^^^ at the moment when T reaches T™^^) and interestingly within the age 
range for which we have many observations available. 

Armed with these considerations, we turn now to compare the predictions 
of the minimal cooling paradigm with data. 
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Figure 25: Neutron star cooling trajectories varying the magnitude of the neu- 
tron ^P-F2 gap. See text for description. 



8.2 Comparison of IVIinimal Cooling With Observations 

Fig. |26| shows ranges of predicted thermal evolutions of a canonical neutron star 
of 1.4Mo buih with the EOS of Akmal & Pandharipande [16|. The density 
dependence of the symmetry energy in this EOS precludes the occurrence of 
the direct Urea process in a 1.4M0 neutron star. Furthermore, this EOS does 
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Figure 26: Comparison of two minimal cooling scenarios with observational 
data. The neutron ■^P-F2 gaps employed are shown in Fig. 12 See text for 
description. Adapted from the results of [52] 
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not permit hyperons, kaon condensates, pions or deconfined quark matter, all 
of which could lead to enhanced cooling. It was shown in [ST] that all EOS's 
compatible with the restrictions imposed by the minimal cooling scenario yield 
almost identical predictions. Moreover, within the range of neutron ^P-F2 pair- 
ing gaps explored, for which is reached at a density that is smaller than 



the central density of any neutron star (see right panel of Fig. 12 1, the neutron 
star mass of any model also has very little effect 

The various possible assumptions about the neutron ^Sq gap have only a 



small effect on the early crust thermal relaxation phase, phase (B) of Fig. 20 
Changes in the proton ^Sq gap produce significantly differing effects, but dom- 
inant effects are due to changes in the neutron ■^P-F2 gaps. The full range of 
possibilities is (almost) covered by considering the two neutron '^P-F2 gaps in 



and "b" of Fig. [12] and is depicted in the two panels of Fig. 26 The 
various curves in the grey shaded areas show the uncertainty in the predictions 
due to lack of knowledge of the envelope compositions, which can range from 
pure heavy elements ("HEE") to pure light elements ("LEE"). In agreement 



with the results of Fig.[25]and the discussion of the previous subsection, the large 
gap "b" results in warmer stars than the smaller gap "a" during the neutrino 
cooling era. 

Comparing model tracks with observations, we can conclude, following [51] 
and [S2|, that if the neutron 3p-F2 gap is of small size {T^^"" - 10^ K) as that of 
model "a" all neutron stars with detected thermal emission are compatible with 
the minimal cooling scenario with the exception of the oldest objects (labelled 



10, 12, and 13 in Fig. 26, see the discussion at the end of Sec. 7.5). These 



older stars are candidates for the presence of internal heating, i.e., the "H" 



term in Eq. ( 50 ) jl23| , a possibility that can easily be incorporated within the 
minimal cooling scenario. In addition, it is found that young neutron stars 
must have heterogeneous envelope compositions: some must have light-element 
compositions and some must have heavy-element compositions, as noted above. 
The updated comparison by j52j more precisely quantifies the required size 
of this gap. In the notations used in Fig. [24] maximal compatibility of the 
minimal scenario requires the neutron ^P-F2 gap to satisfy ^ 5 x 10^ K 

and T™" < 2 x 10* K. The constraint on T™'^'^, which determines when the PBF 
process will be triggered, derives from the necessity of having a PBF enhanced 
already acting in the youngest observed stars at ages ^ 10^ yrs. On the 
other side, the constraint on the low T™'", which determines when core neutrino 
emission will be strongly suppressed, including the PBF process, is obtained 
by the requirement that should not be strongly suppressed before the star 
reaches an age of a few times lO'' yrs. This constraint assures compatibility with 
the low Tg of objects 5 and 6. Considering the smooth kp, and thus density, 
dependence of gaps, the constraint on T™'" likely prevents T'^^^ from being too 
large: gaps with T™^'^ > 10^ K usually also have T™" > 10*^ K. 



The lofFe group, in their version of minimal cooling |136l 1137] , assumed neutron ^P-F2 
gaps that are small at pnuc and grow rapidly at high densities. Consequently, T™**" is reached 
in the center of the star and stars of increasing mass have increasing values of resulting 
in a mass dependence of the cooling evolution. 
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In the case the neutron ^P-F2 gap is as large as our model "b" , one reaches 
the opposite conclusion that minimal cooling cannot explain about half of the 
young isolated stars, implying the occurence of some fast neutrino process. The 
same conclusion is also reached in the case this gap is very small. As seen in 
the left panel of Fig. [25j a vanishing gap leads to cooling trajectories very close 
to those obtained using large gaps during the neutrino cooling era. 

Only one object, the pulsar in the supernova remnant CTA 1 (object "D" in 
the figure) for which only an upper limit on is available, stands out as being 
significantly below all predictions of minimal cooling. 

Finally, irrespective of the magnitude of the neutron ■^P-F2 gap, if any of 
the six upper limits marked "a" to "f in Fig. [26] are, in fact, neutron stars, 
they can only be explained by enhanced cooling. As it is unlikely that all 
of these remnants contain black holes, since the predicted neutron star/black 
hole abundance ratio from gravitational collapse supernovae is not small, these 
objects provide additional evidence in favor of enhanced cooling in some neutron 
starj^ As the nuclear symmetry energy likely increases continuously with 
density, larger-mass stars, which have larger central densities, likely have larger 
proton fractions and a greater probability of enhanced cooling. The dichotomy 
between minimal cooling and enhanced cooling might simply be due to a critical 
neutron star mass above which the direct Urea process can operate before being 
quenched by superfiuidity. 



8.3 Conclusions from the Minimal Cooling Paradigm 

The study of the minimal cooling scenario allows progress beyond the conclu- 
sions reached in Sec. |7.8[ but we are still faced with a clear dichotomy: 

(1) In the case the neutron ■^P-F2 gap is small in size (satisfying T™'" < 
2 X 10'* K and ^ 5 x 10^ K), neutron stars undergoing fast neutrino 

cooling must be relatively rare as among the dozen of known young isolated 
neutron stars, only one candidate, the pulsar in the supernova remnant CTAl, 
exists. This candidate is augmented by possible neutron stars represented by 



one or more of the(green) upper limits in Fig. 26 Nevertheless, the total number 
of fast cooling candidates remains small. 

(2) In the case the neutron ^P-F2 gap is either larger, or smaller than in (1), 
a much larger fraction of neutron stars appear to undergo fast neutrino cooling. 

Both conclusions are extremely interesting, but more information is needed 
to choose the (hopefully) correct one. Such information may be at hand and is 
presented in the next section. These conclusions also have implications in terms 
of neutron star masses. Neutron star mass measurements are only available 
for about 3 dozen pulsars in binary systems, and these may not be indicative 



^•^ Further evidence for enhanced cooUng is provided by two neutron stars, SAX J1808.4- 
3658 and IH 1905-1-00, in transiently accreting low-mass X-ray binaries. In contrast to the 
the six, yet to be detected, candidates "a" to "f , these two stars are known to be neutron 
stars from their characteristics during the accretion phases, but their thermal emission is not 
detectable after accretion stops, implying extremely fast neutrino emission occurring in their 
cores [T381IT391. 



58 



of the overall neutron star population. However, we can appeal to theoretical 



predictions of the neutron star initial mass function, which are shown in Fig. 27 



The predicted distributions are actually similar to the observed distributions 
from binary neutron stars |140j . In the case that conclusion (1) above is correct, 
the critical mass M^^ for fast neutrino cooling would appear to be much larger 
than the average mass, perhaps in the range 1.6 — I.^Mq. Alternatively, if 
conclusion (2) above is correct, ~ IAMq, around the average mass. 
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Figure 27: The initial mass function of neutron stars as predicted by stellar 
evolution models. The continuous line shows results from [141) and the dotted 
line is adapted from [TU]. Figure from [T5] . 



9 Cassiopeia A and its Cooling Neutron Star 

The Cassiopeia A (Cas A) supernova remnant (SNR) was discovered in radio in 
1947 and is the second brightest radio source in the sky (after the Sun). It has 
since then been observed at almost all wavelengths. Very likely, the supernova 
was observed by the first Astronomer Royal, John Flamsteed |142) who, on 
August 16, 1680, when describing the stars in the Cassiopeia constellation, 
listed the star "3 Cassiopeia" at a position almost coincident with the supernova 
remnant. This star had never been reported previously, and was never to be 
seen again, until August 1999 when the first light observation of Chandra found 



a point source in the center of the remnant (see Fig. 28 ) . 

The distance to the SNR is 3. 4^°'^ kpc [T33j, and the direct observation, 
by the Hubble Space Telescope, of the remnant expansion implies a birth in 
the second half of the 17**^ century |144j and supports Flamsteed's observation. 
These observations give a present age of 333 yrs for the neutron star in Cas A. 
The optical spectrum of the supernova has been observed through its light echo 
from scattering of the original light by a cloud of interstellar dust and shows 
the supernova was of type lib ^145^ . The progenitor was thus a red supergiant 
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Figure 28: The Cassiopeia supernova remnant in X-rays: first light of Chcindra, 
August 1999. (The neutron star is highlighted by the authors.) Image from 
©NASA/CXC/SAO. 

that had lost most of its hydrogen envelope, with an estimated zero age main 
sequence (ZAMS) mass of 16 to 20 Mq [1461 EZl EM even up to 25 Mq 
in the case of a binary system |149j . This implies a relatively massive neutron 
star, i.e. likely ^1.4 Mq [149] . The large amount of circumstellar material 
associated with mass loss from its massive progenitor could have diminished 
its visibility from Earth and could explain why it wasn't as bright as the two 
renaissance supernovae, Kepler's SN 1572 and Tycho's SN 1604. 

The soft X-ray spectrum of the point source in the center of the SNR in 
Cas A is thermal, but its interpretation has been challenging [150] . With a 
known distance, a measurement of the temperature implies a measurement of 
the star's radius, but spectral fits with a blackbody or a H atmosphere model 
resulted in an estimated radius of 0.5 and 2 km, respectively. This suggests 
a hot spot, but that should lead to spin-induced variations in the X-ray flux, 
which are not observed. It was only in 2009 that a successful model was found: a 
non-magnetizetp^ C atmosphere, which implies a stellar radius between 8 to 18 
km [151] . Models with heavier elements, or a blackbody, produced significantly 
poorer fits. With the C model, and analyzing 5 observations of the SNR, Heinke 
& Ho |152) found that the inferred neutron star surface temperature had dropped 
by 4% from 2000 to 2009, from 2.12 to 2.04 x 10^ K, and the observed flux had 
decreased by 21%. The neutron star in Cas A, the youngest known neutron 
star, is the flrst one whose cooling has been observed in real time! 

^*There is, to date, no evidence for the presence of a significant magnetic field in the Cas 
A neutron star. 
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9.1 Superfluid Neutrons in the Core of Cas A 



The Chandra observations of Cas A, together with its known distance, imply 
that the photon himinosity of the neutron star is 



10^^ erg s"^ 



(62) 



With a measured Tg°° ~ 2 x lO'^ K [HT], we deduce an internal T ~ 4 x 10* K 
from Eq. (56). The star's total specific heat is thus Cy ~ 4x 10"^* erg (from 

4% [152| gives a change of core 
Assuming 



Fig.[T8|or Eq. (157|)). The observed AT^ /T° 
temperature AT/T ~ 8% over a ten years period since T 



the observed cooling corresponds to the global cooling of the neutron star, its 
thermal energy loss is 



E^i^ = Cvf ~ (4 X 10^** erg K"^) x (0.1 K s" 



4 X 10^^ erg s"^ 



(63) 



which is 3-4 orders of magnitude larger than L-^ ' 



For a young neutron star, 
neutrinos are the prime candidates to induce such a large energy loss. 

The cooling rate of this neutron star is so large that it must be a transitory 
event, which was initiated only recently. Otherwise it would have cooled so much 
it would probably now be unobservable. Something critical occurred recently 
within this star! "Something critical" for a cooling neutron star points toward a 
critical temperature, and a phase transition is a good candidate. The previous 
section highlighted that a phase of accelerated cooling occurs when the neutron 
^P-F2 pairing phase transition is triggered. With T^^^ ~ 5 x 10* K, a transitory 
cooling episode can occur at an age ~ 300 yrs as shown in the right panel of 



Fig. 25 and in Fig. 29 




Figure 29: Similar to the right panel of Fig. [25j but showing the observed age 
and temperature of the Cas A neutron star (the star) and its consistency with 
Tc ~ 5 X 10* K for the magnitude of the neutron '^P-F2 gap. Taken from [55] . 

The interpretation that the observed rapid cooling of the neutron star in 
Cas A was triggered by the recent onset of the neutron •^P-F2 superfluid phase 
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transition and consequent neutrino emission from the formation and breaking 
of Cooper pairs in the neutron superfluid was recently proposed in [S^ and, 
independently, in |153j . 



A simple analytical model 



The simple analytical solution of Eq. ( 58 ) gives some insight into the observed 
behavior. When T > Tc = T^"^, but < Tq, the star follows the asymptotic 
"MU trajectory", Tg = {tmu /^Y^^ , and when T reaches Tc-, at time t = tc-, 



the neutrino luminosity suddenly inc reases (see Fig. 24). Despite the compli- 
cated T dependence of e^^^, Eq. (46), the resulting luminosity, once integrated 



over the entire core (also, aided by the bell shape of the Tc{kp) curve), is well 
approximated by a power law in the T regime in which some thick shell of 
neutrons is going through the phase transition. If we write 

(64) 

gives 



-PBF 



9 > 



with / ^ 10, the solution of Eq. (58), replacing to by tc and Tq by Tc 



Tc 



(65) 



[l + f{t~tc)/tcV/^ ' V 

Thus, at late times, the asymptotic "PBF trajectory" is a factor /^^^^ ~ 0.7 
lower than the "MU trajectory" while the time-scale for the transit from the 
"MU trajectory" to the asymptotic "PBF trajectory", Tg = (tmu//*)^^*^, is 
T'TR = tc/ f ■ This behavior is shown schematically in the left panel of Fig. 30 

Since the initiation of the rapid transit, when the neutron star left the MU 
trajectory, must have occurred recently, i.e. tc — (0.5 — 0.9) x 333 yrs, we obtain 



Ti rrfU 



10^ K {TMu/k 



a/6 



5 X 10*K. 



(66) 

This is the first important result from this simple analytical model: given the 
observation of rapid cooling, the inferred value of Tc depends only on the known 
age, t = 333 yrs, of this neutron star and the value of tmu ^ 1 year. This 
result is also in good agreement with the numerical models of Fig. 29 The 1 /6 



exponent in Eq. ( 66 ) confirms it to be a robust result without much dependence 



on microphysical details. 



9.2 Superconducting Protons in the Core of Cas A 

How "rapid" is the observed "rapid cooling" of the Cas A neutron star can be 
quantified by the observed slope 

rflnTf° 



-^obs 



dint 



1.4 



(67) 



from the six data points reported in |1521 1153) . This slope is much larger than 
the one of the "MU trajectory" , 
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(68) 
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Figure 30: Schematic cooling trajectories (heavy curves) showing the effect 
of superconductivity. Left paneh Without superconductivity, T initially follows 
the modified Urea (MU) trajectory, Tg = (tmu/^)^^^ until T reaches the neutron 
critical temperature Tc at time tc- The pair breaking and formation (PBF) 
process turns on and the neutrino luminosity L,^ abruptly increases by a factor 
/. Thereafter, T rapidly transits, on a time scale r^j^ = tc/ f toward the PBF 
trajectory, Tg — (tmu//0^^^- Empirically, this transition takes a time ^ 2tc- 
Right panel: When protons are superconducting, the initial evolution follows 
a superconducting-suppressed modified Urea (SMU) path. For the transit to 
start at a time t'c ~ tc the trajectory requires T^ > Tc- The early transit has 
a shorter time scale, r^-^ = t'cl f with /' ;» /, and a significantly larger slope. 
The late time evolution converges to that of the left panel. The left (right) panel 



corresponds to models in Fig. ^ (Fig. 31 and 32 1. From [53] 



In contrast, the "transit trajectory" of Eq. (65) gives a much larger slope 



1 f f T Y t 



which has a maximum value //12. Nevertheless, to match the observed slope 
Sobs J / ^ 10 is required. 

The value of / considered above, / = L^^^ /L^^ ~ 10, arose from a freely 
acting modified Urea neutrino emission. However, despite the many theoretical 
uncertainties discussed in Sec. |4j there seems little doubt that proton supercon- 
ductivity exists around a few times Pnuc, as illustrated in Fig. |11[ Moreover, 
expected values of T^'^^ are somewhat larger than 10^ K, implying that protons 
were likely already superconducting in some part of the core of the Cas A neu- 
tron star when neutron anisotropic superfluidity set in. If such was the case, the 
previous neutrino luminosity of this star was due to a proton-pairing-suppressed 
modified Urea process with Lf,^^ < ^f^. This implies a much higher relative 
efficiency of the neutron PBF process, i.e., /' = iP^^/i^MU > / ^ lo. The 
resulting transit from a "SMU trajectory" to the "PBF trajectory" is depicted 
in the right panel of Fig. [30] and exhibits a transit slope enhanced by a factor 
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Figure 31: A fit to observations of the neutron star in Cas A assuming recent 
onset of neutron '^P-F2 superfluidity and PBF cooling. The IAMq model shown 
assumes the APR EOS [16] with a C atmosphere [m]. With model "CCDK" 
of Fig. 11 (T'™^'' ^ 10^ K), proton are superconducting from early times. The 
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K. This neutron superfluid phase transition triggered the PBF process that 
results in a sudden cooling of the neutron star. Observations [1521 1153j suggest 
I d In /dint I ~ 1.4, shown in the inset. Two dotted curves with Tc = and 
1 X 10^ K, respectively, illustrate the sensitivity to Tq- Figure adapted from 



With the above considerations, a very good fit to the observations can be 
obtained, as shown in Fig. 31 implying a maximum neutron '^P-F2 pairing 



Tc = 5.5 X 10® K along with superconducting protons with a larger T™^^. Very 
similar results were independently obtained by |153p^ This observation of 
the cooling of the youngest known neutron star is unique and its interpretation 
potentially imposes very strong constraints on the physics of ultra-dense matter. 

The requirement that protons became superconducting before the onset of 
neutron superfluidity places a constraint not only on the proton ^Sq pairing but 
also on the neutron star mass. If the neutron star mass is not too large, mod- 
els show that proton superconductivity does not extend to very high densities. 



^^Prcsumably, at the time of the onset of proton superconductivity, another PBF episode 
had occurred, but, as the PBF process for singlet pairing is much less efficient than for triplet 
pairing, and also because protons are much less abundant than neutrons, it did not result in 
a significant cooling of the star. Moreover, this cooling occurred on timescales much smaller 
than the crustal thermal timescale. 

These authors, however, assumed that proton superconductivity extends to very high 
densities, with Tcp ^5 (2 — 3) X 10^ K in the whole neutron star core. This results in a very 
strong suppression of prior to the onset of neutron ^P-F2 superfluidity, and a neutron 
star much warmer than in the model of 1531 . from which a larger T^'^ ~ (6 — 9) X 10* K is 
deduced. However, strong proton superconductivity at high densities is, at present time, not 
supported by the microscopic models presented in Fig. 
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Fig. 32 illustrates the sensitivity to the neutron star mass. A better understand- 
ing of the progenitor of Cas A and its expected neutron star mass will prove 
important in validating this scenario. 




Figure 32: Dependence of the slope s = \dlnT^ /dlnt\ of the cooling curve 
on the star mass at t = 330 years: s = 1.4,0.9, and 0.5 for M = 1.3,1.6, and 
1.9Mq, respectively (from [55]). 

The neutron superfluid explanation for the rapid cooling of the neutron star 
in Cas A fits well within the minimal cooling scenario [51 > 52 . As described in 



Sec. 8.2 maximal compatibility of the minimal cooling scenario |52j with data 
required the neutron ^P-F2 gap to have T™""" > 5x 10^. This lower limit on T™'^'' 
was deduced for compatibility with the measured Tg of the youngest neutron 
stars of age ~ 10^ yrs The upper limit on T™'" "^^-^ deduced for compatibility 
with the oldest middle-aged stars as, e.g., the Vela pulsar. The compatibility 
of the neutron ■^P-F2 gap inferred from the cooling data of the neutron star in 
Cas A with observations of other isolated neutron stars is confimed in Fig. [33) 
The only marked difference between Fig. [33] and the left panel of Fig. [26] is the 
occurrence of the rapid cooling phase at ages ~ 300 yrs to ^ 500 yrs, due to the 
reduced value of T^'^^ in the former, 0.55 x 10^ K compared to 1.0 x 10^ K. 

9.3 Alternative Scenarios 

Alternative explanations for the observed rapid cooling of Cas A have been 
proposed. One could consider observed cooling of this star to be due to a sig- 
nificantly longer thermal relaxation timescale in the crust or core than assumed 



in [551 1153j . In such a case, the estimate of Eq. (63), which assumes the star 
to be isothermal, becomes invalid. In all models shown in this chapter, the 
crust thermal relaxation occurs on a timescale of a few decades. However, if 
the crust thermal conductivity is, in fact, significantly smaller, it is conceivable 
that the observed rapid cooling corresonds to the thermal relaxation of the crust 
(see, e.g., [154| ). Such a low crust thermal conductivity is, however, in conflict 
with the observed crust relaxation time in transiently accreting neutron stars 
[1551 1156j and is based on the assumption that the crust is in an amorphous 
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Neutron ^P-Fj gap: "a2" (Xr" =5.5xlO^K) 
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Figure 33: Comparison of the cooling scenario of Fig. [31] with data for isolated 
neutron stars. As in Fig. |26[ the various lines show the effect of varying the 
amount of light elements in the envelope (from |126| ). 

solid state instead of a crystaline one, a possibility that is not supported by 
microscopic studies |157l I158j . 

Similarly, the core thermal relaxation time may be much larger than usually 
considered. For example, Blaschke et al. [159J have proposed that the inner core 
of the star cools rapidly and that it also takes a few hundreds years for the star 
to become isothermal. The latter time is when the rapid decrease of would 
be observed. This scenario requires that the core thermal conductivity be lower 
than usually considered, by a factor 4 or larger, and also requires that neutrons 
do not form a superfluid until the star is much colder. This scenario, based on the 
"Medium-Modified Urea" neutrino emission process [7H I104| , is also compatible 
with the cooling data, but only if the suppression of core conductivity is adjusted 
to fit the observed cooling of the neutron star in Cas A. 

More work is required to confront these alternative possibilities with other 
facets of neutron star phenomenology. 

10 Dynamical Evolution of Neutron Stars 

Pulsars are rotating neutron stars whose spin rates are generally observed to 
be decreasing. The regularity of pulsars is outstanding; as timekeepers they 
rival atomic clocks. Although the pulses are remarkably regular, the time be- 
tween pulses slowly but predictably increases. Their spin-downs are attributed 
to magnetic dipole radiation [1601 1161) - the conversion of rotational energy 
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into electromagnetic energy - and the pulsar wind from ejection of the mag- 
netospheric plasma |162j . That the observed evolution of pulsar spins exhibits 
evidence for both core and crustal superfluidity is the subject of this section. 



10.1 Magnetic Dipole Pulsar Model 

Dimensional analysis gives us a simple estimate of the pulsar energy losses. 
Energy and angular momentum are irreversibly lost when either the magnetic 
field or the plasma reaches the light cylinder where co-rotation with the pulsar 
imphes a speed equal to c and whose radius is thus Ric = c/fl ~ 50(P/1 msec) 
km, n being the spin rate and P — 2tt/Q the spin period. The magnetic field at 
the light cylinder is Bjc = {R/Ric)^Bs assuming a dipolar field with a strength 
Bs at the stellar surface. Just writing that an energy density Bf^/Air is lost at 
the speed of light c from a sphere of area AnRf^, one obtains 



EpsHO,\{R^B,,fn\ (70) 

Energy loss by in vacuum magneto-dipolar radiation gives the same result with 
just an extra factor (2/3)sin^Q;, a being the angle between the rotationnal 
axis and the dipolar moment. An aligned rotator, i.e. with a = 0, will not 
spin-down from magneto-dipolar radiation but rather by plasma ejection [162] . 
Numerical and consistent calculations of the energy loss from the ejected magne- 
tospheric plasma have been possible only recently ( '163' and |164j for a review) 



and the result is that the total EpsR is given by Eq. ( 70 ) with an extra factor 
~ (1 + sin^ a). Considering the rotational energy of a uniformly rotating sphere 
Erot = /ri^/2, where / is the star's moment of inertia, the pulsar's spin-down 
is determined by equating Ej-ot with —EpsR, giving 

ti = -jn'^ (71) 

with K ~ R^B^/c^, assuming that / and Bs remain constant. This also provides 
an estimate of B^, probably reliable within a factor of a few, 



1/2 

n^iR^ ' ^ ' ' 0.01 s V T, 



3.2 X 10-(PP)V2 G . 1.5 X 10-^ C-^) G, 



(72) 

we assumed AI ~ 1.4 M0 and i? ~ 12 km. It is traditional to also deduce an 
observable characteristic pulsar age 



and an observable braking index 



Tc = -^ (73) 
20 ^ ' 



n = ^ (74) 
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In the case of the magnetic dipole model and assuming both / and Bs are 
constant, n — 3. 

The ATFN cataloguej^ [165] hsts 1759 pulsars (as of December 2012) with 
measured values of P. The measurement of a second derivative Cl, allowing 
the determination of n, requires a very accurate long term and smooth fit of 
Cl and is much more difficult due in part to timing noise and to glitches. To 
date, there are only 11 published values of n [166], all with values less than 
3 and ranging from —1.5 to +2.91 (see Table 10.2). Anomalous values n < 3 



may be due to a growing magnetic field |167] and/or to a growing fraction of the 
star's core becoming superfluid |168) . The latter could be due to the onset and 
growth of ^So p and '^P-F2 n components (further evidence for this stems from 
observations of the cooling of the neutron star in the Cassiopeia A supernova 
remnant, which is discussed in Sec. |9]). 



10.2 Anomalous Braking Indices 

Consider the evolution of the spin of a pulsar if we allow changes in the surface 
field Bg and the moment of inertia /. Assume the star has both normal matter 
and superfluid matter with moments of inertia /„ and Ig , respectively. Although 
the onset of superfluidity is unlikely to significantly change the total moment 
of inertia, the portion of the star that is superfluid may be considered to spin 
frictionlessly and its spin rate Qs to remain constant. The observed spin rate of 
the star is that of the star's surface, which is composed of normal matter. We 
therefore have 



I^In+Is, 



nr. 



(75) 



With these two components, the rate of change of angular momentum in the 
magnetic dipole model is 



dt 



[I A + is^s] = inii + in {n - = -Kn\ 



Taking a time derivative of the above, one finds 



n = 3 + 4tc 



In Bs 
In Bs 



- 1 



In 
In 



(76) 



(77) 



Since fi^ — ^ fi, we can drop the last term in Eq. (77 1. Thus, the growth 



of either the surface magnetic field and/or the superfluid component leads to 
n < 3. In fact, the growth of a core superfluid component could lead to field 
ejection from the core and an increase in the surface field strength. Note that 
there has been a long history of study of the opposite possibility of decay of 
pulsar magnetic fields. However, statistical studies of the pulsar population 
|169j no longer support this idea and the present consensus is that magnetic 



On-line version: http://www.atnf . csiro . au/research/pulsar/psrcat/ 
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J1734-3333 


1.1693 


8.13 


52.2 


0.9 






B 1737-30 


0.6069 


20.6 


17 




0.79 


1.2 


J1747-2958 


0.0988 


25.5 


2.49 


< 1.3 






B1757-24 (G5.4-1.2) 


0.1249 


15.5 


4.04 


-1 


1.35 


1.5 


B1758-23 (W28?) 


0.4158 


58.3 


6.93 




0.24 


1.0 


B1800-21 (G8.7-0.1(?)) 


0.1367 


15.8 


4.29 


2 


1.57 


1.8 


B1823-13 


0.1015 


21.4 


2.8 


2 


0.78 


1.2 


J1846-0258 (Kes75) 


0.3266 


0.728 


48.8 


2.6 






B1930+22 


0.1445 


39.8 


2.92 




0.95 


2.8 


B2229+6114 (G106.6+2.9) 


0.0516 


10.5 


2.03 




0.63 


0.5 



Table 2: Observed pulsar properties ("SNR" is the associated supernova rem- 
nant). Values of spin period P, characteristic age Tc [Eq. ( 73 )] and dipole surface 
field strength [Eq. ( [72|)] r ounded from the ATFN catalogue [BS]. Values of 
the braking index n [Eq. (74)] rounded from |166j . Values of A [Eq. (79)] taken 
from 



field decay is significant only in the case of super-strong fields as in magnetars 
or after a long phase of accretion in a binary system. 

Nevertheless, observed values n < 3 do not prove that superfluidity exists. 
The growth of the surface field of a pulsar can be obtained in the case that the 
initial magnetic field was partially buried, e.g., by post-supernova hypercritical 
accretion |170l 11711 1172| . into the stellar interior and is now slowly diffusing 
outwards [167] . Detailed numerical models show that observed values of n can be 
reproduced |173j without invoking superfluidity. Alternatively, magnetospheric 
currents may deform the field into a non-dipolar geometry |174l I175j that results 
in n < 3 (but n <C 3 is unreachable in these models and n < is impossible). 
The truth may be a superposition of these three types of arguments: /„ < 0, 
Bg > 0, and non-dipolar fields. 



10.3 Glitches 

The second indication of superfluidity in neutron stars stems from the fact that 
many pulsars exhibit sporadic spin jumps, or glitches (see Fig. 34). The Jodrell 
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Figure 34: Schematic illustration of a Vela pulsar glitch. The observed pulsar 
spin frequency v has been corrected by the average accumulated spindown since 
1969.0 {vq = ii{t — 1969.0)) to make the glitch easier to discern. The glitch 
is observed as a relatively sudden jump Ai/ in spin frequency. The frequency 
afterwards relaxes roughly exponentially on a timescale t to a new spindown 
line (dotted) which has a higher average frequency {\ — Q)/S.v than the pre-glitch 
spindown line v = (dashed). 



Bank glitch catalogue]^ hsts 420 glitches (as of December 2012) [T7Z]- These 
are thought to represent angular momentum transfer between the crust and a 
liquid, possibly superfluid, interior |178|ll79[ll80j . As the star's crust spins down 
under the influence of magnetic torque, differential rotation develops between 
the crust (and whatever other parts of the star are tightly coupled to it) and a 
different portion of the interior containing a superfluid (recall in the above that 
J7s = 0), but possibly a superfluid component different from those in the core. 
The now more rapidly rotating (superfluid) component then acts as an angular 
momentum reservoir which exerts a spin-up torque on the crust as a result of 
an instability. The Vela pulsar, one of the most active glitching pulsars, glitches 
about every 3 years, with a fractional change in the rotation rate averaging a 



part in a million [1811 1182| , as shown in Fig. 35 

The stocastic nature of glitches implies that they represent a self-regulating 
instability for which the star prepares over a waiting interval. The amount of 
angular momentum transferred is observed to increase linearly with time, as 
shown in Fig.[35]for Vela and for the X-ray pulsar J0537-6910. Assume that the 
moment of inertia of the superfluid component responsible for glitches is Ig and 
that of the crust and whatever parts of the star are tightly coupled to it is Ic- 
In an individual glitch i, an amount of angular momentum JcAfii is transferred 



On-line version: http://ww.jb.man.ac.uk/pulsar/glitches.html 
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Figure 35: The accumulated relative angular velocity 10^^ J^i Af2i/f2 as a func- 
tion of time in days for the X-ray pulsar J0537-6910 and for Vela (B0833-45). 
Arii is the change in angular velocity of a single glitch. Linear fits are indi- 
cated, with respective slopes [ Eq. A ^ 2A x 10"^ d"^ and 1.91 x 10~^ 
d-\ (From [ m\ ) 



to the crust, where Afti is the observed jump in angular velocity. Over the 
interval ti between the last glitch and the current glitch, the star spins down by 
the amount AQ = —nti and the total angular momentum differentially stored 
in the superfluid component is IgAQ. These two angular momenta are equal, 
and over a total observed time tots = J2i have 

Ig > 2t,AIc, (78) 

where the observed quantity A is related to the accumulated jumps, 

A= {^An,/^t-,\. (79) 

In the case of the Vela pulsar, the magnitude of A implies that Ig/Ic = ^.TcA ^ 
0.016 [183]. Minimum values for this ratio for other pulsars are given in the last 
column of Table [TOl 

Since the moment of inertia in the inner crust (with the n ^Sq type superfluid) 
is Isc and satisfies 

Isc ^ 0.04/, (80) 

where / is the star's total moment of inertia, glitches can be naturally explained 
by the inner-crust superfluid since Ig ^ Isc- Although the magnitudes of in- 
dividual glitches varies somewhat, the maximum-sized glitches have stable and 
limited sizes. This would be difficult to explain if glitches originated from the 
inner core's superfluid component, whose associated moment of inertia is a much 
larger fraction of the star's total. 

Recently it has been shown that most of the neutron mass of the dripped 
neutrons in the inner crust is entrained by Bragg scattering with the nuclear 
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lattice, effectively increasing the neutron mass by factors of 4-5 |1841I178] . With 
this entrainment, Eq. (78) becomes, for Vela, 



Ig > 0.08/e, (81) 

precluding an inner crust superfluid explanation for Ig if ~ /. Entrainment 
results in most of the neutron fluid spinning down with the crust, and the 
unentrained conduction neutrons cannot accumulate angular momentum at a 
high-enough rate to produce the largest observed glitches |185j . 

B. Link jl85j has pointed out, however, that Ic is likely to be smaller than 
/, perhaps considerably so, if the outer core decouples from the crust over 
timescales of weeks to years. In fact, this is the observed timescale r of the 



post-glitch relaxation (see Fig. 34 1 , which may simply represent the dynamical 
recovery of the outer core. For older pulsars, which are cooler, the relaxation 
timescale is long, of the order of years, and this could explain the nearly step-like 
behavior of many pulsar glitches. As a consequence, for example, if /c — //2, 



Eq. (81) and (80) regain consistency. The upshot is that spin glitches could 
then originate in either the inner crust or the outer core. The observed limited 
maximum magnitude of glitches retains its natural explanation in terms of the 
small and regular sizes of neutron star crusts as opposed to the wider variations 
in sizes of core superfluid regions. 

It should be noted that the above is at odds with the results of Sidery and 
Alpar |186| . who obtain relaxation times of only about 2 days. This discrepancy 
is the subject of ongoing discussions. 



11 Discussion and Conclusions 

In this chapter, the influence of pairing, leading to neutron superfluidity and 
proton (and quark) superconductivity, on key observables of neutron stars is de- 
scribed. The observables include aspects of their thermal evolution, composed 
of surface temperatures, cooling rates, and ages, and dynamical evolution, com- 
prised of pulsar spin-down characteristics and glitch information. The major 
effects of pairing on the thermal evolution of an isolated neutron star are the 
quenching of neutrino emissions and the reduction of specific heat of the paired 
fermions, be they nucleons, hyperons or deconfined quarks, in their cores. How- 
ever, the onset of pairing also triggers short episodes of increased neutrino emis- 
sion through pair breaking and formation (PBF) processes when the ambient 
temperature falls below the superfluid critical temperature. The PBF emission 
of neutrino pairs through weak interactions of strongly interacting particles is 
unique to dense neutron-star matter, as a similar phenomenon does not occur 
in nuclei for ^Sq paired fermions as it is forbidden on the basis of symmetry. The 
major effects of pairing on dynamical evolution include the reduction of the so- 
called braking index, a measure of spin deceleration, of pulsars, and the possible 
triggering of glitches from weak coupling of superfluid vortices to neutron star 
crusts. 
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The sizes and density dependences of superfluid and superconducting gaps 
play a crucial role. Pairing is observed in nuclei as an energy difference between 
even-even and odd-even nuclei with a typical magnitude, ranging from 0.5 to 3 
MeV, that decreases with atomic number. The basic cause of pairing in nuclei 
is due to the attractive interaction between neutrons in the spin S = channel. 
Gaps of similar magnitude for nucleon pairing in neutron stars are expected. 
Both spin-singlet -'^Sq (at lower densities) and spin-triplet (at higher den- 
sities) configurations appear possible from scattering phase shifts arguments. 
Pairing appears as a gap A in the single particle energy spectrum, leading to 
a strong suppression (^ e~^^/^) of both specific heat and neutrino cmissivi- 
ties at low temperatures. In the simple BCS approximation, however, the gap 
depends exponentialy on the pairing potential. Hence, uncertainties associated 
with in-medium effects of strong-interactions at high density have prevented a 
consensus about the sizes and density dependences of gaps from being reached. 
In addition, the nature and abundance of possible candidates for pairing (nu- 
clcons and strangeness-bearing hyperons, quarks, Rose condensates) are also 
uncertain. 

In the neutron star crust, the density of unbound neutrons is high enough 

that -'^Sq pairing is expected to occur. Most theoretical models suggest that the 
associated gap disappears at neutron densities higher than that of the core-crust 
interface, so it is confined to the crust. The latest calculations indicate a max- 
imum gap magnitude of about 1.5 x 10^" K. Since proton-neutron correlations 
reduce the effective mass of the proton below that of the neutron, the size of 
the proton ^Sq gap is smaller than that of the neutron. "Unbound" protons 
exist only at densities greater than pnuc where nuclei disappear, so proton su- 
perconductivity in the spin-singlet state is expected to exist from the core-crust 
boundary to deep into the core once temperatures fall below a few times 10^ 
K. The "'^So neutron gap vanishes close to the core-crust interface and the dom- 
inant pairing for neutrons in the core occurs in the anisotropic ^P-F2 channel. 
Uncertainties in the size and density range for this gap are larger than for ^Sq 
gaps, however, with maximum magnitudes ranging from a few times lO'^ K to 
a few times 10^ K. 

The greatest influence of pairing will be on the thermal evolution of neutron 

stars. The occurrence of pairing leads to three important effects for neutron 
star cooling: alteration and eventual suppression of nucleon specific heats, sup- 
pression of neutrino emissivities, and triggering of PBF neutrino emission for 
temperatures just below the critical temperature. Presently, 13 isolated neu- 
tron stars with thermal spectra in the soft X-ray band have been identified. In 
addition, there are four pulsars with detected X-rays but with only upper lim- 
its to thermal emission. Finally, there arc six gravitational-collapse supernova 
remnants that might contain neutron stars (if not, then black holes), but no 
detected thermal emission as yet. Atmospheric modeling of the thermal sources 
yields estimates of surface temperatures, and together with age estimates, al- 
lows these stars to be compared to theoretical cooling models. Most observed 
sources are younger than a million years, during which time they cool primarily 
through neutrino emission. 
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Neutrino cooling can be either very fast (i.e., enhanced), or relatively slow. 
Enhanced cooling occurs by way of the direct Urea process on nucleons, hyper- 
ons, Bose condensates, or deconfined quark matter. It is allowed when energy 
and momentum can be conserved with 3 or fewer degenerate fermions involved. 
If additional "bystander" nucleons are required to conserve momentum, neu- 
trino emission is suppressed by about a factor of a million. If pairing is not 
present and enhanced neutrino emission does not occur, it is found that several 
observed neutron stars are too cold to match cooling models. Therefore, an 
important first conclusion is that either pairing or enhanced neutrino emission 
must occur in some neutron stars. On the other hand, if pairing is not present 
and enhanced neutrino emission does occur, all observed neutron stars are too 
hot to match cooling models. Thus, a second conclusion is that in the presence 
of enhanced cooling, superfluidity and/or superconductivity must occur. How- 
ever, as there are many combinations of pairing gap sizes and extents, neutron 
star masses and envelope compositions, and enhanced cooling candidates that 
can match observations, it is not possible to determine either gap properties or 
the specific enhanced cooling reactions involved. Therefore, for example, neu- 
tron star cooling cannot lead to an unambiguous detection of hyperons, Bose 
condensates or deconfined quark matter in the interior of neutron stars at this 
time. 

A final scenario, known as the minimal cooling paradigm, assumes enhanced 
neutrino emission does not occur, but allows pairing. In this case, all observed 
stars, with the possible excption of the pulsar in the supernova remnant CTA 1, 
for which only an upper temperature limit is available, and the six undetected 
neutron stars in gravitational-collapse supernova remnants, none of which is 
certain to be a neutron star at all, can be fit by theoretical cooling models. 
However, fitting observations requires two constraints on neutron stars: 

• The neutron '^P-F2 gap must have a maximum critical temperature T^'^ 
larger than ^ 0.5 x 10^ K, but unlikely larger than 1.5 — 2 x 10^ K; and 

• some, but not all, neutron stars must have envelopes composed of light 
elements (H/He/C) and some, but not all, must have envelopes composed 
of heavy elements. 

The assumed mass of a neutron star has only a minor influence on the results 
of the minimal cooling scenario. The source CTA 1 and some of the undetected 
neutron stars in gravitational-collapse supernova remnants might then be candi- 
dates for enhanced neutrino emission via the direct Urea process. One obvious 
way in which this could occur is that, since there is a density threshold for the 
onset of the direct Urea process, only neutron stars above a critical mass could 
participate. From observed and theoretically-predicted neutron star mass dis- 
tributions, this critical mass is estimated to be in the range 1.6 — l.SM©. This 
result is itself supported by recent experimental and observational restrictions 
on how fast the symmetry energy of nuclear matter can increase with density 
|140j . which suggest that the nucleon direct Urea threshold density is much 
greater than pnuc- 
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The recently detected rapid cooling of the neutron star in Cas A provides an 
unprecedented opportunity to overcome the dilemma of deciding whether "exot- 
ica" are needed, or if the minimal cooling paradigm is sufficient, to account for 
observational cooling data. The rapidity of the stellar cooling in Cas A remnant 
points to both the onset of core neutron superfluidity within the last few decades 
and the prior existence of core proton superconductivity with a larger critical 
temperature. The PBF process, with its time-dependent burst of neutrinos at 
the critical temperature, is central to the success of model calculations fitting 
both Cas A cooling observations and the entire body of data for other observed 
isolated neutron stars. Furthermore, the Cas A cooling observations can be fit 
only if T^'^ ~ 5 — 6 x 10* K, precisely in the range independently found for 
the minimal cooling scenario based on other isolated neutron star cooling data. 
Sherlock Holmes would not have deemed this a coincidence. "Enhanced" pro- 
cesses appear to be ruled out in the case of Cas A, as its star would presently 
be too cold. Further observations of the cooling of the neutron star in Cas A 
can confirm these conclusions. 

Final indications of the existence of pairing in neutron stars is afforded by ob- 
servations of the deceleration of the spin frequency and glitches of pulsars. The 
standard paradigm to explain the pulsar mechanism is based on the magnetic 
rotating dipole model, for which it can be shown that the spin-down parameter 
n = — nO/n^ = 3. All non-accreting pulsars for which n has been determined 
have 71 < 3, which can be understood either if the surface dipolar magnetic field 
is increasing in strength or if there is a growing superluid or superconducting 
component within the neutron star's core. An increasing surface field could, 
in fact, be due to the expulsion of magnetic flux from the core by supercon- 
ductors. The leading model for pulsar glitches is that a more-rapidly rotating, 
and essentially frictionless, superfluid in the crust sporadically transfers angular 
momentum to the more slowly-rotating crust and the parts of the star strongly 
coupled to it. 

In conclusion, there is abundant evidence that the ultimate "high-temper- 
ature" superfluid (or superconductor) exists in nearly every neutron star. Un- 
mistakably, superfluidity and superconductivity exhibited by the constituents 
of neutron stars provide avenues by which the observed cooling behaviors of 
these stars can be explained. Theoretical predictions of pairing gaps being as 
yet uncertain, astronomical data are pointing to a path forward in their deter- 
minations. A remarkable recent development is that, despite the nearly three- 
order-of-magnitude theoretical uncertainty in the size of the neutron ^P-F2 gap, 
and the factor of 3 uncertainty in the size of the proton ^Sq gap, observations of 
the neutron star in Cas A now appear to restrict their magnitudes to remark- 
ably small ranges. It will be interesting to see if these considerable restrictions 
translate into a more complete understanding of the pairing interaction between 
nucleons and provide insights into other aspects of condensed matter and nuclear 
physics. Continued X-ray observations of this star and other isolated sources, 
and, hopefully, the discovery of additional cooling neutron stars, will further 
enhance these efforts. 
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